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Some turbulent solutions of the unaveraged Navier-Stokes equations (equa- 
tions of fluid motion) are reviewed. Those equations are solved numerically 
in order to study the nonlinear physics of incompressible turbulent flow. 
Initial three-dimensional cosine velocity fluctuations and periodic boundary 
conditions are used in most of the work considered. The three components of 
the mean-square velocity fluctuations are initially equal for the conditions 
chosen. The resulting solutions show characteristics of turbulence, such as 
the linear and nonlinear excitation of small-scale fluctuations. For the 
stronger fluctuations the initially nonrandom flow develops into an apparently 
random turbulence. Thus randomness or turbulence can arise as a consequence 
of the structure of the Navier-Stokes equations. The cases considered include 
turbulence that is statistically homogeneous or inhomogeneous and isotropic or 
anisotropic. A mean shear is present in some cases. A statistically steady- 
state turbulence is obtained by using a spatially periodic body force. Various 
turbulence processes, including the transfer of energy between eddy sizes and 
between directional components and the production, dissipation, and spatial 
diffusion of turbulence, are considered. It is concluded that the physical 
processes occurring in turbulence can be profitably studied numerically. 


I. INTRODUCTION 

Nearly all of the flows occurring in nature, as well as those that are 
manmade, are turbulent. For instance, the boundary between a column of rising 
smoke and the surrounding atmosphere is generally irregular and contains a 
range of scales of motion, indicating the presence of turbulence. The atmos- 
phere itself is usually turbulent, as shown by the irregular appearance of many 
of the clouds present in it. Jets, wakes, astrophysical flows, and flows over 
surfaces are commonly turbulent, as is the region downstream of a grid in a 
wind tunnel or downstream of a waterfall. In general, turbulent flows are the 
rule and laminar flows the exception. 

Because of the importance and challenge of the turbulence problem a great 
deal of research has been done over the past century. Basic ideas have been 
set forth, for instance, in papers by Reynolds (1883, 1895), Taylor (1921, 
1935), von Karman (1937a, 1937b), and Heisenberg (1948). That work, together 
with more recent research, is discussed in books by Batchelor (1953), Hinze 
(1975), Frost and Moulden (1977), and others. 

In spite of considerable research activity there is no general deductive 
theory of high-Reynolds-number (strong) turbulence. (Reynolds number is 
defined as the product of a velocity and a length divided by the kinematic 
viscosity of the fluid. It is a measure of the ratio of inertial to viscous 
effects.) Most of the analytical theories depend on a closure assumption for 
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a hierarchy of averaged equations. 1 This Immediately calls Into question the 
appropriateness of referring to the analytical theories as deductive, except 
at low Reynolds numbers or In the final stage of decay, where nonlinear effects 
are small (Batchelor, 1953; Delssler, 1977). 

One way In which a closure assumption can be avoided Is by closure by 
specification of sufficient random Initial conditions (Delssler, 1979). That 
method can successfully predict turbulence decay, and In the sense that the 
evolution of all Initially specified quantities can be calculated, gives a 
complete solution. To use It, however, the Initial conditions must be fully 
turbulent, and a large amount of Initial data Is required to satisfactorily 
specify the Initial turbulence. The method does not seem capable of extension 
to some cases of turbulence maintained by mean gradients, where the effect of 
Initial conditions may eventually become negligible (e.g., In fully developed 
turbulent pipe flow) . 

In view of the foregoing comments It seems desirable to consider numerical 
solutions of the unaveraged Navler-Stokes equations that display features of 
turbulence. Numerical methods and computers can be considered as tools for the 
solution of equations, just as can Fourier transforms and series expansions. 

It might be pointed out that It Is more appropriate to refer to a numerical 
solution of the unaveraged equations as deductive than It Is to so refer to 
most of the analytical theories, which are based on averaged equations and 
require closure assumptions. Moreover, most of the analytical theories are so 
complicated that a large amount of numerical work Is required to obtain results 
from them. Attempts to obtain analytical solutions of the unaveraged equations 
have not been successful, mostly because of the nonlinearity of those equa- 
tions. Herring (1973) mentions that the simplest turbulence theory Is just the 
Navler-Stokes equations. Since most turbulence calculations are numerical 
anyway, no Insight Is lost by considering direct Integration of the Navler- 
Stokes equations forward In time, starting from some suitable Initial data. 

Numerical solution (or numerical simulation) has sometimes been called 
experiment. It seems, at least to this writer, that there Is an Important 
difference between numerical solution and experiment as generally practiced. 

The former uses directly, and attempts to solve, a given set of constitutive 
equations, In this case the Navler-Stokes equations. The latter ordinarily 
does not, although both methods may arrive at the same result If the constitu- 
tive equations are congruous with the portion of nature to which they are 
applied. In general, It appears that experiment works directly, with nature, 
whereas numerical solution works with a set of constitutive equations that 
should represent at least a portion of nature. 

Several numerical solutions of the unaveraged equations have appeared that 
use a spectrum of random Initial fluctuations (e.g., Orszag and Patterson, 

1972; Clark et al . , 1979; Rogallo, 1981; and Felerelsen et al . , 1982). These 
studies, which appear to demonstrate the feasibility of carrying out turbulent 
solutions with present-day computing equipment, represent major advances. 


Hhe hierarchy of correlation (averaged) equations obtained from the 
unaveraged Navler-Stokes equations Is unclosed because of the nonlinearity of 
the latter. That Is, there are more unknowns than equations, so that a 
closure assumption Is required to obtain a solution. 
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Because of the difficulty of specifying realistic turbulent Initial con- 
ditions (experimentally or analytically), It may be more appropriate to Ini- 
tially specify a simple regular fluctuation with a single length scale (as 
actually occurs downstream of a grid In a wind tunnel). This should be better 
for studying the development of small-scale fluctuations (and of turbulence In 
general) than would a spectrum of Initial fluctuations, since for the latter, 
small-scale fluctuations are already present In the Initial flow. Moreover, 
much higher Reynolds number flows can be calculated with a given numerical grid 
when a single length scale Is Initially present, at least for early and moder- 
ate times. Taylor and Green (1937) and others (e.g., Delssler, 1970a; 

Van Dyke, 1975; Corrsln and Kollman, 1977; and Delssler and Rosenbaum, 1973) 
have used a perturbation series to calculate the nonlinear development of 
higher harmonics from lower ones, but the calculations could not be carried 
very far In time. In these analyses the directional components of the Initial 
fluctuation Intensity were not equal. Orszag and Fateman (Orszag, 1977a) have 
recently used Taylor and Green's Initial conditions and obtained a numerical 
solution for higher Reynolds numbers and longer times. The Invlscld (Infinite 
Reynolds number) case was Investigated In some detail by Betchov and Szewczyk 
(1978) . 

The present review considers the nonlinear physics of turbulence numeri- 
cally. Although the Initial conditions used herein are nonrandom, the flow at 
higher Reynolds numbers breaks up into an apparently random turbulence. Unlike 
the problem of Taylor and Green all three of the directional components of the 
mean-square velocity fluctuations are equal at the Initial time. In the 
absence of mean shear they are also equal at later times. Taylor and Green's 
directional components, on the other hand, do not approach equality, even at 
long times (Orszag, 1977a). 

To study the processes In turbulence, first some background on the basic 
fluid flow and turbulence equations Is given In section II and on numerical 
methods and solutions In section III. Four cases of turbulence are then con- 
sidered, starting In section IV with the simplest one, In which mean gradients 
are absent (Delssler, 1981a). In this case no energy sources are present and 
the turbulence decays freely. (By contrast the presence of mean gradients 
would Imply energy sources In the flow.) Here (In section IV) one can study 
viscous dissipation and the nonlinear transfer of energy between wave numbers 
or eddy sizes, as well as the randomization of the flow. Next, In section V a 
uniform mean shear Is applied to study turbulence production and maintenance, 
and the linear and nonlinear transfer of energy between wave numbers and 
between directional components (Delssler, 1981c). The transfer of energy 
between wave numbers (both linear and nonlinear) Is manifested by the creation 
of small-scale structure In the turbulence. Then, In section VI, the spatial 
diffusion of the Inhomogeneous turbulence In a developing shear layer Is con- 
sidered (Delssler, 1982). Finally, by using a spatially periodic body force a 
turbulence that Is statistically steady state at long times Is studied In sec- 
tion VII. The first three of these cases have also been studied, but with a 
spectrum of random Initial fluctuations and In some cases with an assumption 
for the small eddies. In Orszag and Patterson (1972), Clark et al. (1979), 
Rogallo (1981), Shannan et al. (1975), and Cain et al. (1981). Here we will 
confine ourselves to the development of turbulence from nonrandom Initial 
conditions with a single length scale. 

One of the problems In the numerical study of turbulence Is that of accu- 
racy, because of the small scale of some of the turbulent eddies. As the 
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Reynolds number (strength) of the turbulence or the time Increases, smaller 
eddies are generated. No matter how small the numerical mesh size, one can 
always pick a Reynolds number or time large enough that the results will be 
quantitatively Inaccurate. One way of Improving the accuracy Is by extrapola- 
tion to zero mesh size, as will be done here. The effectiveness of that pro- 
cedure depends to some extent on the accuracy of the unextrapolated solution. 

If the solution has to be extrapolated too far, the results may not be accu- 
rate. Another popular method (not used here) Is to model eddies smaller than 
the grid spacing (subgrid modeling) (e.g., Smagorlnsky, et al., 1963; 

Deardorff, 1970; Clark et al., 1979; and Ferzlger, 1977). This method requires 
an empirical Input, although not as great a one as that for full modeling of 
the averaged equations. One might think of subgrid modeling as a useful crutch 
that can be phased out as numerical resolution Improves. However, when It Is 
used. It Is sometimes difficult to tell which effects come from the equations 
of motion and which result from the subgrid modeling. Slggla (1981) has 
recently considered the converse problem; he made a numerical study of the 
small-scale eddies In which he modeled the larger ones. 

Here we are mainly concerned with physical processes and trends, rather 
than with highly accurate numerical results (possibly unattainable at very high 
Reynolds numbers). Of course some degree of accuracy Is necessary; otherwise 
we will not even be able to calculate trends. As the numerical mesh size 
decreases, quantitative differences In the results might be obtained. It Is 
to be hoped, however, that the results will not be qualitatively different. 
Results to date Indicate that to be the case. 2 

Other relevant review articles are given by Orszag (1977b), Schumann, et 
al. (1980), Eckman (1981), and Ott (1981). 


II. BASIC EQUATIONS AND CONCEPTS 
A. The Unaveraged Equations 

Turbulent flows of a great many liquids and gases obey the Navler-Stokes 
equations. Those equations assume that the fluid Is Newtonian (stress propor- 
tional to strain rate) and that It can be considered a continuum. The latter 


2 The attainment of accurate quantitative results appears to be a ques- 
tion of Improvement of computers and of numerical methods. If state-of-the-art 
numerical methods and computers are used, good quantitative as well as quali- 
tative results can already be obtained, at least for low and moderate Reynolds 
numbers. Orszag and Patera (1981) (as well as Moln and Kim (1982) using sub- 
grid modeling) made significant numerical calculations of the velocity profile 
In the wall region of fully developed turbulent channel flow. The results 
agreed reasonably well with experiment, showing a wall transition region and a 
fully turbulent region In which the velocity varies as the logarithm of dis- 
tance from the wall. The advent of high-speed computers and efficient numeri- 
cal algorithms may be making possible for the first time the use of the 
Navler-Stokes equations In the solution of a wide range of realistic 
(turbulent) fluid-flow problems. 
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Is usually a good assumption because In most cases Intermolecular lengths are 
much smaller than the smallest significant turbulent eddies. 

1 . Equations In terms of Instantaneous quantities 

The Navler-Stokes and continuity equations for constant fluid properties 
(Including Incompressibility) 3 * 4 can be written as (e.g., Batchelor, 1967; 
or Delssler, 1976) 



at 
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The subscripts can take on the values 1, 2, and 3, and a repeated subscript In 
a term Indicates a summation, with the subscript successively taking on the 

values 1, 2, and 3. The quantity ui Is an Instantaneous velocity compo- 
nent, x-| Is a space coordinate, t Is the time, p Is the density, v the 
kinematic viscosity, and p Is the Instantaneous pressure. Equations (1) and 
(2) are, respectively, statements of the conservation of momentum and of mass. 
In order to obtain an explicit equation for the pressure, we take the diver- 
gence of equation (1) and apply the continuity equation (2) to get 
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In the remainder of the paper It will usually be convenient to use 
equations (1) and (3) rather than (1) and (2). Equations (1) (1 = 1, 2, 3) 
and (3) constitute a set of four equations In the four unknowns u^ and p. 
Since they are for Instantaneous velocities and pressures, they should apply 
to turbulent as well as to laminar flows, subject to the restrictions mentioned 
at the beginning of this section. The Navler-Stokes equations have been known 
for more than a century, but their use In turbulent flows, other than In a 
schematic sense, has been restricted by a lack of ability to obtain solutions. 
Now, with advances In computers and numerical methods the situation appears 
somewhat brighter. 


3 The continuity equation is sometimes Included In the Navler-Stokes 
equations . 

4 Most turbulence studies have been carried out for constant properties, 
for simplicity. That flow Is realistic If the turbulence velocities are rea- 
sonably low compared with the velocity of sound, and If temperature gradients 
are not large. 
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The fundamental turbulence problem Is an Initial-value problem. That is, 
given initial values for the u^ as functions of position, a value for v, 
and suitable boundary conditions, equations (1) and (3) should be sufficient 
for calculating the u^ and p/p as functions of time and position. The 
initial and boundary conditions used herein are specified in section III. 

To Interpret the terms in equation (1), it is convenient to multiply it 
through by p and by the stationary volume element dx-j dx 2 dX 3 . Then the 
term on the left side of the equation is the time rate of change of momentum in 
the element /ou-j dx-| dx 2 dx 3 - This rate of change is contributed to by the 
terms on the right side of the equation. The first term on the right side, a 
nonlinear inertia term, Is the net rate of flow of momentum into the element 
through its faces. The next term, also nonlinear, is a pressure-force term 
and gives the net force acting on the element by virtue of the pressure gradi- 
ent in the x^ direction. It is nonlinear because of the nonlinear source 
term on the right side of the Poisson equation for the pressure (eq. (3)). 
Finally, the last term in equation (1), a linear viscous-force term, gives the 
net force acting on the element in the x^ direction by virtue of the 
viscosity. 

2. Equations in terms of mean and fluctuating components 

Following Reynolds (1895) one can break the Instantaneous velocities and 
pressure in equations (1) to (3) into mean and fluctuating (or turbulent) 
components; that is, set 


U 1 » U 1 


+ u. 


(4) 


and 


p = p + p 


(5) 


where 


u 1 = p = 0 (6) 

U 1 = ^ (7) 

and 

P - P (8) 

The overbars designate averaged values. ^ Equation (2) becomes, on using 
equations (4), (6), and (7), 




(9) 
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which shows that both the fluctuating and mean velocity components obey con- 
tinuity. Equations (1) and (3) become, on using equations (4) to (9), taking 
averages, and subtracting the averaged equations from the unaveraged ones, 
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Equations (10) and (11) will be used to study the processes in turbulence, but 
not for computational purposes (except In linearized cases). The first four 
terms of equation (10) and the first two of equation (11) look like the terms 
In equations (1) and (3), although their meanings are exactly the same only 
If = P = 0 (eqs. (4) and (5)). The first three terms on the right side of 
equation (10), which contribute to au^/at, can still be Interpreted as an 
Inertia-force (or turbulence self-interaction) term, a pressure-force term, and 
a viscous-force term. The remaining terms are, respectively, a turbulence 
production term, a mean-flow convection term, and a mean-turbulent-stress term, 
which appears when the turbulence Is statistically Inhomogeneous (when mean 


turbulence quantities such as u^u^ are functions of position). (The reasons 

calling the production and convection terms as such will perhaps become clearer 
when the equivalent terms In the averaged equations (eqs. (14) and (15)) are 
discussed.) It will be seen that when the mean velocity gradient Is not zero, 
the term -U^ au^/ax^ generates a small-scale structure In the turbulence by 
vortex stretching. The nonlinear self-interaction term -a(u^U| < )/ax| c also 
produces a small-scale structure and In addition produces randomization of the 
flow. These effects will be considered In sections IV to VI. The Poisson 
equation for the pressure fluctuation (eq. (11)) has three source terms: a 

nonlinear term, a mean-gradient term, and a mean-turbulent-stress term, which 
appears when the turbulence Is Inhomogeneous. 


We have defined an Inhomogeneous turbulence as one In which averaged 
turbulence quantities are functions of position. Thus, a homogeneous turbu- 
lence Is one In which averaged turbulence quantities are not functions of 
position. For Instance, In homogenous turbulence 


5 For the most general flows the average Is usually an ensemble average 
over a large number of macroscoplcally Identical flows ( 1 . e . , mean quantities, 
but not fluctuating quantities, are the same In all the flows). In most cases, 
however, statistical uniformity or statlonarity with respect to one or more 
coordinates, or with respect to time, obtains. Then, the average Is taken with 
respect to the one or more coordinates or with respect to time. According to 
the ergodlc theorem those averages are the same as the ensemble average If the 
flow Is turbulent. In section V (uniform mean shear) three-dimensional spatial 
averages are used, even when the periodic boundary conditions Introduce some 
local Inhomogeneity Into the fluctuations. Those averages still have meaning 
since their values are Independent of the position of the boundaries of the 
cycle over which the averages are taken. 
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Similar statements apply to other averaged turbulence quantities in a 
homogenous turbulence. 

B. Averaged Equations 

Although the averaged equations will not be solved numerically because 
they do not form a closed set, they are useful for studying the processes In 
turbulence. 


1 . Equations for mean flow 


First consider the equations obtained by averaging 
equations (1) and (3) and using equations (4) to (9 ). 5 
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These equations look like equations (1) and (3) with Instantaneous values 
replaced by average values, but with the Important difference that an extra 

term Involving the quantity now appears In each of the equations. 

These terms arise from the nonlinear velocity terms in equations (1) and (3) 
and are a manifestation of the closure problem of turbulencel. If those 
terms were absent, equations (12) and (13) could be solved, and turbulent flows 
would be no more difficult to calculate than laminar ones. Note that terms In 
equations (12) and (13) that contain lower-case letters (other than x's) are 
turbulent terms. 


The form of equation (12) suggests that the quantity -pu u augments the 
viscous stress p«3Ui/3X|<. Since It involves the fluctuatlngor turbulent 
velocity components u^ and u|<, we Interpret It as a turbulent or 

Reynolds stress. For Instance -pu u will. In the presence of a mean-velocity 
gradient 3U-|/3X2, act like a shear stress on an X 1 -X 3 plane. In the 
presence of 3 Ui/ 3 X 2 , U] will more likely be negative than positive when U 2 
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Is positive, so that will have a nonzero negative value. The quantity 

-pu,u_ can be compared with a viscous shear stress obtained In the kinetic 
theory of gases, where u-j and U 2 are now molecular, rather than macros Dpi c 

2 

velocity components. Similarly pu, will act like a normal stress on an 
X 2 -X 3 plane (similar to a normal stress or pressure obtained In klnet.c 
theory, where u-j Is again taken to be a molecular velocity component). 


2. One-point correlation equations 
a. Construction of equations 


We can construct equations for the undetermined quantities u^u^ In 

equations (12) and (13) from the evolution equation for u^ (eq. (10)) and 
a similar equation for the component uj: 
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Multiply equation (10) by uj and the preceding equation for Uj by 
u^ , add the two equations, and average. ^ This gives, using continuity 
(eq. (9)), 
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Setting 1 = j and using continuity, we get, for the rate of change of the 
kinetic energy per unit mass, 
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As In equations (12) and (13), quantities In equations (14) and (15) that 
contain lower-case letters (other than x's) are turbulent quantities. 
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The one-point correlation equation (eq. (14)) gives an expression for the 


rate of change of u^u^ that might be used In conjunction with equations (12) 
and (13). But the situation with respect to closure Is now worse than It was 
before. Whereas without equation (14) we only had to determine u^u^ , with It 
we have to determine quantities like u^u^u^, pu^, p 3u^/ax^, and 
(au^/ax^MaUj/ax^) . We might use equation (11) to obtain the pressure corre- 
lations, but that would only Introduce more unknowns. However, equations (14) 
and (15) are useful for studying the processes In turbulence in that most of 
the terms have clear physical meanings. Moreover, we will be able to calcu- 
late terms In those equations from our numerical solutions. 

b. Physical Interpretation of terms 


As In the case of equation (1) It Is helpful, for purposes of Interpre- 
tation, to multiply the terms In equations (14) and (15) through by p and 
by a volume element dx dx ? dx . Then the term on the left side of 


equation (14) or (15) gives the time rate of change of 


pu iY 


or of the 


kinetic energy pu,u,/2, within the element. This rate of change Is contrib- 
uted to by the termson the right side of the equations. The first of those 
terms Is equal to the net work done on the element by turbulent stresses act- 
ing In conjunction with mean-velocity gradients. It Is therefore called a tur- 


bulence production term; It equals the rate of production of pu^ or of 

pu,u,/2 within the volume element by work done on the element. A somewhat 
abbreviated Interpretation suggested by the form of the term, which Is often 


given, Is that It represents work done on the turbulent stress pu.u, by the 
mean-velocity gradient. 1 


The next term In each of the equations describes the convection or net 
flow of turbulence or turbulent energy Into a volume element by the mean veloc- 
ity Ufc. It moves the turbulence bodily, rather than doing work on It by 
deforming It, as In the case of the production term. It vanishes when either 
U|< Is zero (no mean flow) or when the turbulence Is homogenous 

(u.u. ji u.u, (x.)). In the latter case there Is no accumulation of turbulence 
1 j 1 j k 

within a volume element, even with a mean flow. 


The next three terms In equation (14) and In equation (15) also vanish 
for homogeneous turbulence. Since they do not contain the mean velocity, they 
do not convect or move the turbulence bodily. Therefore we Interpret them as 
diffusion terms, which diffuse net turbulence from one part of the turbulent 
field to another by virtue of Its Inhomogeneity. The pressure-veloclty- 
gradlent terms In equation (14) drop out of the contracted equation (15) 
because of continuity (eq. (9)). Therefore, they contribute nothing to the 


rate of change of the total energy u^u^/2, but they can distribute the energy 

2 

energy among the three directional components u^/2 ( no sum on 1)- The last 

term In equations (14) and (15) Is the viscous dissipation term, which dis- 
sipates turbulence by the presence of fluctuating velocity gradients. 
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3. Two-point correlation and spectral equations 


To consider the transfer of turbulence between eddy sizes or wave numbers 
(spectral transfer), we must obtain two-point correlation equations. Terms 
related to that process do not appear In the single-point equations ((14) and 
(15)). To obtain two-point equations, we use In addition to equations (10) 
and (11), written at the point P, the following equations written at the 
point P : 
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Multiplying equation (10) by u' and equation (16) by u^, adding, taking 

averages, 5 and using equation (16) and the fact that quantities at one point 
are Independent of the position of the other point result In 
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To simplify the equations, the turbulence will be considered homogeneous 
(correlations Independent of x^). (See Hlnze (1975) or Delssler (1961) 

for Inhomogeneous equations.) Then the two-point correlations (e.g., u.u' ) 
will , be functions only of r, = x! - x,, so that a/ax, = -a/ar, and ^ 
a/ax^ = a/ar^ . In addition, homogeneity requires that au^/ax^ be a con- 
stant, so that au,/ax, = au,/ax, and u k - u k = r, au,/ax,. (If au,/ax, 
were a function or x^, we Joula not remove arl x'j dependency from 1 J 
eq. (18).) Equation (18) becomes 
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The equations for the pressure-velocity correlations obtained from 
equations (11) and (17) are, for homogeneous turbulence. 
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a. Interpretation of terms 

The terms on the right side of the Poisson equations (eqs. (20) and (21)) 
are source terms associated with the mean velocity and with triple correlations 
arising from the nonlinear velocity terms In equations (11) and (17). The 
right side of the two-point equation (eq. (19)) contains, as In the case of the 
single-point equation (eq. (14)), turbulence production, directional distribu- 
tion (pressure-velocity), and viscous terms. There are no diffusion terms In 
equation (19) since the turbulence Is homogeneous. The last two terms In that 
equation are new terms that do not have counterparts In the single-point equa- 
tions. To Interpret them, we convert equation (19) to spectral form by taking 
Its Fourier transform. Thus, define the following three-dimensional Fourier 
transforms: 


Vj (r) = J «p^j( K)e 


it 


(22) 


_3 

3r 


u^ 1 


^ isj ■ £ 


(23) 


dll. 3 U.U! 

—k r JL.3 

3x„ 8. 3r. 




f T ij( V K - r d^ 


(24) 


and 


3 

3r. 


u i u i u k - u A u j ■ / V“ )e ' K ' r 


(25) 


where <p^, n^, T^, and T^ are, respectively, Fourier transforms of the 

quantities on the left sides of thq defining equations (eqs. (22) to (25)), 
k Is a wave-number vector, and dK = d*, d* 2 dic 3 . Physical Interpretations 
of the Fourier transforms defined by equations (22) to (25) follow. 

Equation (19) becomes, on taking Its Fourier transform, 


3U, 


3U 


->2 


3t V 1J ' " v kJ 3X. " “Mk 3X. + p n 1j (,c) “ 2uk ’ij + T 1 j + T lj 


(26) 


If we let r = 0 In equation (22), we have 
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( 21 ) 


u 1 u j = J m V 1J dK 


so that we can Interpret <p., as a spectral component of u u,. As In 
equation (19), terms on theVlght side of equation (26) can be J Interpreted as 
production, directional distribution, and viscous dissipation terms. They 

contribute to the rate of change of a spectral component <p^ of u^. 

To Interpret the term T^ In equation (26), we let r = 0 In 
equation (24). This gives 



(28) 


Thus, T' gives zero total contribution to the rate of change of u^. But 
It can distribute, or transfer between wave numbers or eddy sizes, spectral 
contributions to u^u^ . So T^, which Is proportional to aU^/ax^, Is 

Interpreted as a mean-gradient spectral transfer term (Delssler, 1961). The 
term -r^au^uj/ar^ aU^/ax^ In equation (19) Is therefore the Fourier trans- 
form of a mean-gradient spectral transfer term. 

To Interpret the last term In equation (26), we use equation (25), where 
we note that a/ar = a/ax' = -a/ax , that quantities at one point are Inde- 
pendent of the position of k the other point, and that continuity (eq. (9)) 
holds. Equation (25) becomes 



1*.r 


dK 



u 1 u j u k 



u 1 u k u 


j 


au! au u. 
= -Vk ax£ - u j 


If we let r = 0, equation (29) becomes 




u iYk - 0 


(29) 


(30) 


since we have assumed homogeneity of the turbulence. Thus, as In the case of 

It can, how- 


T 1 T 

i r u 


contributes nothing to the rate of change of u^. 


ever, transfer spectral components of u.u. from one part of wave-number 
space to another. So we Interpret T,, n as a spectral transfer term asso- 
ciated with turbulence self-lnteractlin (as opposed to Interaction between 
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turbulence and mean gradients). The term -(a/ar^Hu^u^ - u^u^) In 

equation (19) Is therefore the Fourier transform of a self-interaction spectral 
transfer term. Although the turbulence has been assumed to be homogeneous for 
Interpreting and as transfer terms, It has been shown that similar 

Interpretations apply when J the turbulence Is Inhomogeneous (Delssler, 1981b). & 

The transfer of turbulence from one part of wave-number space to another, 
or from one eddy size to another, produces a wide range of scales of motion In 
most turbulent flows as Illustrated by the numerical solutions In sections IV 
and V (e.g., figs. 2(c), 8, and 24). The state of affairs Is neatly summarized 
In a nonmathematlcal way by a poem written long before equations (26), (28), 
or (30) was known (Richardson, 1922): 


Big whirls have little whirls, 

Which feed on their velocity; 

And little whirls have lesser whirls, 
And so on to viscosity. 


4. Vortlclty and dissipation 

For homogeneous turbulence one can obtain a relation between the viscous 
dissipation term In equation (15) or (19) (1 =^j ) and the vortlclty or swlrl^ln 
the turbulence. The dimensionless vortlclty o> Is defined as the curl of u: 

w(x) = Vxu(x) (31 ) 


or 


"l = G 1jk ax, 


where Is the alternating tensor. Then 


(32) 


w k“k = c 


au^ au ft au^ au^ au^ au^ 

1 jk e Pmk ax^ ax^ (6 1P 6 jm " 6 1m 6 jP ) ax^ i)T “ 


au^ au 


ax 


j ax 


J 


ax^ ax^ 


(33) 


But, because 



(34) 


6 The fact that terms related to T., or Ti, do not appear In the one- 
point eq. (14) for Inhomogeneous turbulence Indicates that T^ and T^ do 

not contribute to auTuT/at. Thus, even for Inhomogeneous turbulence they can 
only transfer turbulence from one part of wave-number space to another. 
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(35) 


for homogeneous turbulence. 


au^ au^ 

3Xj ax^ 


2 - 
a u 

i!i 

ax^ 

3x j ‘ 

(33) 

then 

au 1 

au 1 


(0, G>. = (0 = ^ 

k k ax. ax. 


j j 


£ 

v 


(36) 


Thus, for homogeneous turbulence the mean-square vortlclty Is just the rate of 
viscous dissipation e of turbulent energy divided by the kinematic viscosity 
(eq. (15)). So the more Intense the swirl In the turbulence, the faster It 
dissipates . 


From a closure standpoint we are somewhat better off with the two-point 
equations (eqs. (19) to (21)) than we were with the single-point equation 

(eq. 


(14)), since we no longer have to model terms like 


p aUj/ax^ 


and 


( au^/axjj^ au^/ax^) . However, we still have to evaluate triple-correlation terms 

unless the turbulence Is very weak or unless the Interaction between the turbu- 
lence and the mean flow Is large compared with the turbulence self-interaction. 


5. Remarks 


We shall not discuss here the many schemes that have been proposed for 
closing the averaged equations considered In this section. Instead, we will 
avoid the closure problem by obtaining numerical solutions of the unaveraged 
equations (eqs. (1) and (3)). The Importance of the averaged equations Is 
enhanced by these numerical solutions; by using the solutions of the unaveraged 
equations, terms In the averaged equations that represent various physical 
processes In the turbulence can be calculated. Thus, the averaged equations 
appear to be necessary, or at the least very convenient, for the physical 
Interpretation of the numerical results. 

III. NUMERICAL SOLUTIONS AND METHODS 

A. Initial conditions 

For most of the numerical solutions considered here, the Initial velocity 
fluctuation Is assumed to be given by 

3 

E n -»n 

cos q .x (37) 

n=l 

Then, from equation (4) 


^ e 1jk = 0 when 1, j, and k are not all different. When the sub- 
scripts are all different, e-jjk = +1 when they are In cyclic order and -1 
when they are In acyclic order. 
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(38) 


u 


1 


= a" cos q n .x + 
n=l 


The quantity a" Is an Initial velocity amplitude or Fourier coefficient of 
the velocity fluctuation, q"* n Is an Initial wave-number vector, and Is 
an Initial mean-velocity component. To satisfy the continuity conditions 
(eqs. (2) and (9)), we set 


a"q" = 0 


(39) 


For the present work let 

a] = k(2, ±1, 1), a* = k(l, ±2, 1), a^ = k(l, ±1, 2) 

q] - (-1, ±1, l)/x 0 , q^ = (1, +1, l)/x 0 , q^ = (1, ±1, -l)/x Q 


(40) 


where k has the dimensions of a velocity and determines the Intensity of the 
Initial velocity fluctuation. The quantity xq Is the length scale of the 
Initial velocity fluctuation. The quantities k and xq, together with the 
kinematic viscosity v and equation (40) then determine the Initial Reynolds 

-rl/2 

number Uq Xq/v, since the square of equation (37), averaged over a period, 
~2 

gives Uq. In addition to satisfying the continuity equation (eq. (39)), 
equations (37) and (40) give 



at the Initial time. 8 Thus, equations (37) or (38) and (40) give a particu- 
larly simple Initial condition In that we need specify only one component of 
the mean-square velocity fluctuation. Moreover, for no mean shear they give 
an Isotropic turbulence at later times, as will be seen. Note that It Is nec- 
essary to have at least three terms In the summation In equation (37) or (38) 
to satisfy equation (41). We do not specify an Initial condition for the pres- 
sure because It Is determined by equation (3) and the Initial velocities. 

B. Numerical grid and boundary conditions 

To carry out numerical solutions subject to the Initial condition given by 
equations (37) or (38) and (40), we use a stationary cubical grid with a maxl- 

3 

mum of 32 points and with faces at x^ = x^/Xg = 0 and 2ir. For boundary 
conditions we assume periodicity for the fluctuating quantities. That Is, let 
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bj Is a variable length. 


for b* « 2ir, where b* = bj/Xg, x* = x^/x^, and 
By using equations (4) and (5) these become 

(u.) = (u.) + (U ) - (U ) 

xf=2ir+bt x*=b? xt=2ir+bt x?=bt 

j j j j j j j J 

and 


p x*=2ir+b* = P x*=b* + P x*=2ir+b* P x*=b* 
j j jj j J 

In the present work we assume also that P, given by equation (13), Is 
periodic, so that 


( 44 ) 


(45) 


x?=2ir+b} 
J j 


v*_K* 

X j' D j 


(46) 


and equation (45) becomes 


S*- 2 ’ tb ! ' (4 

These equations are used to calculate numerical derivatives at the boundaries 
of the computational grid. 

C . Numerical solutions 


In carrying out the numerical solutions, we have a choice of solving a 
set of equations containing u,, p, U., and P (eqs. (10) to (13)), or one 
containing u^ and p (eqs. (1) and (3)). The latter set, which Is much 

simpler, Is generally used here. 9 That Is, we solve equations (1) and (3) 
subject to Initial condition (38) and boundary conditions (44) and (47). 

The spatial- and time-differencing schemes (which numerically conserve 
momentum and energy) are essentially those used by Clark et al . (1979). For 
the spatial derivatives In equations (1) and (3) we use centered fourth-order 
difference expressions (e.g., McCormick and Salvadore, 1964). For Instance, 
the fourth-order difference expression used for au^/3x^ Is 


8 The first three terms of eq. (41) apply at all times when there are no 
mean gradients In the flow. 


Ill 



n 



12 k 


(u.) - 8(u.) + 8(u. ) - (u.) 

. n-2 1 n-1 1 n+1 1 n+2j 


where Ax|< Is the grid-point spacing, and the subscripts n, n + 1, etc. 
refer to grid points In the X|< direction. Fourth-order difference expres- 
sions are often considered more efficient than the usual second-order ones 
(Orszag and Israeli, 1974). (Spectral methods devised by Orszag and associates 
are still more efficient but may be somewhat trickier to use.) Centered 
expressions (same number of points on both sides of n, see preceding expres- 
sion) can be used both at Interior grid points and at the boundaries of the 
grid; when n refers to a point on a boundary, values for u. outside of the 

grid, which are required for calculating the numerical derivatives at the 

boundary, are obtained from boundary condition (44). 

For time differencing we use a predictor-corrector method with a second- 
order (leapfrog) predictor and a third-order ( Adams-Moulton) corrector 
(Ceschlno and Kuntzman, 1966). If m represents a time step, and ( ) the 

right side of equation (1^), then at each grid point In space the second-Srder 

leapfrog predictor for u^ at time step m + 1 Is 


(M (1) = (u.) + 2 At(R.) 

m+1 m-1 1 m 


and the third-order Adams-Moulton corrector Is 



( 2 ) 

m+1 


<v + 

m 


At 

12 



+ 8(R ) - (R.) 

m m-1 


where At is the time Increment. The quantity (R^)^ In the preceding 

~ m m+1 

corrector Is calculated by using (u.) v ' In the right side of equation (1), 

~ ( 1 ) m+ ^ 

where (u^) v Is calculated from the leapfrog predictor. Note that the leap- 
frog metho<? + l so-called because It leaps over the time step m) , although unsta- 
ble for all At when used by Itself for Navler-Stokes equations, Is stable 
when used as a predictor. 


The Poisson equation for the pressure (eq. (3)) Is solved directly by a 
fast Fourier transform method. This method of solution was found to preserve 
continuity quite well (v.u « 0) except near the ends of some of the runs for 
no mean gradients, where the solutions began to deteriorate. (Another Indica- 
tion of Incipient solution deterioration near the ends of some of the runs for 
no mean gradients was that equation (41) was no longer accurately satisfied. 8 ) 


^Dividing the velocities and pressures Into mean and fluctuating com- 
ponents Is evidently not advantageous from a computational standpoint except, 
possibly, In linearized cases. 
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Two known types of numerical Instabilities can occur in the present solu- 
tions: a viscous Instability connected with the first and last terms In 

equation (1), which occurs If v At/(Ax k )2 is too large; and a convective 
Instability connected with the first and second terms (or the first and third 
terms through eq. (3)), which occurs If u^ At/ Ax k Is too large. In these 
criteria At, Ax k , and ui are, respectively, a time step, distance step, 
and velocity. Thus, a particular solution should be numerically stable If, for 
a given Ax k , the time step Is sufficiently small. Numerical stability Is 
typically obtained when the solution varies smoothly from time step to time 
step, with no significant breaks In the slope from one step to the next. This 
Is the case for all of the results given here. 

For the present solution very good temporal resolution Is obtained auto- 
matically when At Is sufficiently small to give numerical stability. That 
temporal resolution Is much better than the three-dimensional spatial resolu- 
tion, which Is more severely limited by the storage and power of the computer. 
However, as will be seen (fig. 19), sufficient spatial resolution Is obtained 
to give reasonably accurate averaged results for times not excessively large. 
Some of the averaged results are extrapolated to zero spatial mesh size In an 
effort to obtain more accuracy. The fourth-order method of extrapolation 
(consistent with the fourth-order differencing used here) Is given In Delssler 
(1981a, 1981c). 

IV. HOMOGENEOUS FLUCTUATIONS AND TURBULENCE, NO MEAN FLOW 


For this case (U 
and (3) without the ~ 
numerically 


= 3 u.u /3x. = 0), equations (10) and 
' s overinstantaneous quantities. We 


(11) reduce to (1) 
thus solve 


and 


au i _ a (y k ) iM _ + 

3t ~ 3X k p 3X^ V a * k X k 


(48) 


1 a 2 p 
' 8x a 8x a 


8 <Vk> 
8 Vk 


(49) 


subject to Initial conditions (37) and (40) and boundary conditions (42) and 
(43). In equation (40) for the coefficients In the Initial conditions, we 
choose the first set of signs. 

A. Dimensionless form of equations 

For carrying out the numerical solutions and presenting the results In as 
general and compact a form as possible, we nondlmenslonal 1 ze equations (48) 
and (49) as 


a < u K> 32* t 

3t* = “ 3X* 3X* + 3X*3X* 


(48a) 
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and 


where 


a 2 p* 

3 x j 3x t * 


ax* ax k * 


(49a) 



and xq Is the Initial fluctuation length, which first appeared In 
equation (40). Note that all of the quantities have been nondlmenslonallzed 
by xq and the kinematic viscosity i>. 


The Initial and boundary conditions In equations (37), (40), (42), and 
(43) become In dimensionless form 


u* 


3 

£ 

n=1 


-»n* 

cos q .x 


(37a) 


a]* = k*(2,±l ,1 ) , a 2 * = k*(l,±2,l), a^* = k*(l,±l,2) 

q] * = (-1,±1,1), q^* = (l.+l.l), q^* = (1,+1,-D 

(uf) = (uf) 

x*=2ir+b* x*=b* 


and 


P x 5 =2,r+b J " Px r b 5 

where, In addition to the dimensionless quantities defined above, 


(40a) 


(42a) 


(43a) 


A a " 

> a r 


_n* 
q 1 = 


x o q v 


and 


k* = — k. 

V 


Note that k* has the form of a Reynolds number, where xq Is the length and 
k Is the velocity. The Initial Reynolds number appearing In the figures of 

nV 2 (~i\ n 

this paper, (U qI Xq/v, where (ugj Is an Initial root-mean-square velocity 
component, Is obtained by choosing a value for k* In equations (40a) and 
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space averaging 



and 2, where u,/ 
equations ( 48a j 


2 

u^ over a period, by using equation (37a). The value of 

/v so obtained Is also used In the ordinates of figures 1 

,-^ 1/2 

ui 1 , and u* Is calculated from 

"ting with Initial values from eq. (37a)). 


~ 2\ 1/2 

U Q = U 1 ; 
no {49a) (sta 


B. Development of random fluctuations 


Figure 1 shows the calculated evolution of velocity fluctuations (normal- 
ized by the Initial root-mean- square velocity) at two fixed points In space 
for the Initial Reynolds number shown In the figure (no correction for discre- 
tization error). Since there Is no Input of energy, the fluctuating motion 
eventually decays to a state of rest. In spite of the nonrandom Initial con- 
dition (eq. (37) or (37a)), the velocity fluctuations have the appearance of 
those for a random turbulence. It Is Important to point out that the fluctua- 
tions are not due to numerical Instability since a large number of time steps 
(typically about 20) lie between changes of sign of du^/dt. In connection 
with the high Reynolds number In figure 1 It might be pointed out that Betchov 
and Szewczyk (1978) obtained reasonable turbulent-like numerical results even 
for Infinite Reynolds number for times not excessively large (see also 
Schumann et al., 1980). 

1 . Randomness as sensitivity to Initial conditions 

The dashed curves of u,/(u 0 ] are for Initial conditions perturbed 

approximately 0.1 percent. (The coefficients a^ and q^ given by eq. (40) 

are changed by 0.1 percent of their values.) The perturbed curves follow the 
unperturbed ones for a short time and then depart sharply. Thus, a very small 
perturbation of Initial conditions causes a large change In the values of u^ 
(except near t = 0) . On the other hand, the root-mean-square values of the 
velocities (spatially averaged) decrease smoothly and are unaffected by the 
perturbation of the Initial conditions. All of these features are character- 
istic of turbulence. (The observed sensitivity of the Instantaneous flow to 
small changes In Initial conditions may have unfavorable Implications for 
detailed long-term weather predictions (Lorenz, 1963).) 


We note that the spatially averaged values In figure 1 follow 

approximately the decay law u 2 ~ t -n , where n ~ 2.5. This lies between the 
value for n of 3.3 observed for turbulence downstream of a waterfall (Ling 
and Saad, 1977), and the value of 1.2 generally observed for turbulence gener- 
ated by flow through a grid In a wind tunnel (Uberol, 1963). The decay law Is 
evidently very much dependent on the Initial condition for the turbulence. 


2. Effect of numerical mesh size on randomness 

To get an Idea of the effect of mesh size In the numerical grid on the 

/~ 2\ 1/2 

apparent randomness of the velocity fluctuations, values of u^/IUqI 

at the center of the grid are plotted against t* In figure 2 for three mesh 
sizes. All three of the curves have a random appearance. However, as the 
number of mesh points Increases (as the mesh becomes finer), smaller scale 
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fluctuations are resolved and the randomness appears to be greater. This 
trend Indicates that the observed randomness Is not due to the use of too 
coarse a grid. 

C. Further evidence for randomness and Indications of Isotropic turbulence 

For turbulence the correlation between velocities at two times 

u.j ( t)u.| ( t Q ) should go to zero as the separation of the times t and tg 

Increases. Figure 3 shows that this occurs for the present hlgh-Reynolds- 
number calculations. For true turbulence the correlation should probably 
decrease smoothly with time. This Is nearly the case for the larger tg 
(fig. 3(b)). (tg Is the time In the correlation coefficient In fig. 3 that 
remains fixed as the other (variable) time t Increases.) At early times 
there Is probably some nonrandom structure In the turbulence caused by the 
nonrandom Initial conditions (fig. 3(a)). 

As a further Indication that the hlgh-Reynolds-number flow breaks up 

Into turbulence, we calculate the evolution of the cross correlation u^u,,. 

2 2 2 8 

Although u^ = u 2 = u g at all times, the Initial u^u 2 given by equations 

(37) and (40) Is not zero. However, figure 4 shows that because of the 

apparent randomization of the flow u^u 2 goes to zero as time Increases. 

The fluctuations In the curve at early times (as also In the curve of 
fig. 3(a)) are probably caused by nonrandom structure In the flow at early 


2 

Figures 1 to 4, together with the fact that the three components u^ 

are equal, show that at later times we appear to get a reasonable approxima- 
tion to Isotropic turbulence, although the Initial conditions are nonrandom. 

One of the consequences of Isotropy Is that the cross correlations, say u u , 
are zero, as In figure 4 at later times. These calculations differ from 
others where turbulence was obtained by using random Initial conditions. 

Because the Initial conditions for the present calculations are nonrandom, 
the turbulence must arise as a result of the structure of the Navler-Stokes 
equations. In particular the nonlinear terms play a crucial role. It Is easy 
to verify that If nonlinear terms are neglected for the present case, 
equations (48), (49), (37), and (40) give 

* 

u, = (u.) e" 3t (50) 

1 1 0 

So If the equations are linear, the flow given by the nonrandom Initial 
condition (37) remains nonrandom. The nonlinear terms must be present In 
equations (48) and (49) If the development of turbulence Is to take place. 


10 A11 averaged values would be expected to vary smoothly only for 
highly random fluctuations. 
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Our calculations at lower Reynolds numbers give results that are less like 
turbulence. Thus, the fluctuations develop apparent randomness only at the 
higher Reynolds numbers. 

Calculated values of velocity-derivative skewness factor for a high 
Reynolds number (fig. 5) also appear to be of reasonable magnitude when com- 
pared with experimental values for Isotropic turbulence. Values of that quan- 
tity of the order of 0.4 have been obtained experimentally and have long been 
considered an Indication of true moderately strong or strong Isotropic turbu- 
lence. The fall of f near the end of the curve may be due to a more laminar 
flow there, where fluctuation levels are lower. A longer period of agreement 
with experimental values Is obtained by Clark et al. (1979), apparently because 
they used turbulent Initial conditions, so that the turbulence was already 
partially developed at t = 0. 

D. Origin of the randomness (strange behavior) 

A question remains as to how the nonlinear terms In equations (48) and 
(49) produce the randomness observed In figures 1 to 4. Until recently, It was 
generally assumed that randomness In a turbulent flow Is due to randomness In 
the Initial conditions, to random external fluctuations, or to the presence of 
so many eddies or harmonic components (or of so many degrees of freedom) that 
the Identity of the Individual eddies Is lost (Monin, 1978; Rabinovich, 1978). 
In the present results the first of these Is absent. Concerning the second, 
roundoff errors might be considered a form of external fluctuations. However, 
when the calculations were repeated using double precision, so that roundoff 
errors were reduced by a factor of about 10 6 , the mean-square velocities were 
practically unchanged. The Instantaneous velocities were different, although 
still as random as before. Thus, the effect of a large decrease In roundoff 
errors Is similar to the effect of a small perturbation of the Initial condi- 
tions (fig. 1). Since roundoff errors do not affect the turbulence level or 
the randomness, they cannot be considered a major sustaining cause of the tur- 
bulence or randomness observed here, although they may In some cases affect the 
Initial transition. In the present case the transition Is so rapid that the 
effect appears to be small. 

This leaves only the proliferation of eddies or harmonic components as a 
source of apparent randomness. That might well produce the randomness observed 
In figure 1, since the nonlinear production of harmonics tends to be explosive, 
particularly at high Reynolds numbers (each harmonic component Interacts with 
every other one) (Delssler, 1970a). However, the randomness may be produced, 
at least partially, by strange attractors or, more properly, by analogous 
strange behavior (Eckmann, 1981; Ott, 1981). (We talk about analogous strange 
behavior here, rather than strange attractors, since, strictly speaking, 
strange attractors exist only for steady-state turbulence (section VI I— A) . 

Here, analogous strange behavior refers mainly to apparent randomness In flows 
where a large number of degrees of freedom or harmonic components are not a 
necessary Ingredient and randomness occurs by a loss of hydrodynamic stabil- 
ity.) Lorenz (1963) and others (Monin, 1978; Rabinovich, 1978; Ruelle, 1976; 
Lanford, 1982) have shown that a system of nonlinear ordinary differential 
equations similar to the spatially differenced form of the Navler-Stokes equa- 
tions used here, or to the spectral form of those equations (Orszag, 1977b), 
can develop an apparently random behavior In time as a result of the loss of 
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stability of the solutions.^ In conjunction with this, regions appear In 
the phase space of the system to which solutions are attracted. Randomness 
arises In those regions, which are known as strange attractors, from a haphaz- 
ard movement of the phase point among the neighborhoods of various critical 
points In the phase space (steady-state points of unstable equilibrium where 
au-j/at = 0). (The presence of critical (or fixed) points Is not always con- 
sidered a necessary Ingredient of strange attractors, but randomness or sensi- 
tivity to Initial conditions Is essential. However, the spatially differenced 
and spectral forms of the Navler-Stokes equations do appear to have critical 
points, as do the Lorenz equations.) Unlike randomization by proliferation of 
harmonic components, randomization by strange behavior can occur either with a 
few or with many degrees of freedom. In the present case both processes may be 
Important. At any rate, as mentioned earlier, the results show that the struc- 
ture of the Navler-Stokes equations Is such that apparently random or turbulent 
solutions can arise from nonrandom Initial conditions. With the results from 
the low-order models, In which apparent randomness appears with as few as three 
degrees of freedom (e.g., In the Lorenz equations), the turbulence observed to 
be manufactured by the Navler-Stokes equations should perhaps not come as a 
surprise. 

The presence of strange behavior may be fortunate from a numerical stand- 
point In that It should enable turbulent solutions that are qualitatively cor- 
rect (at least Insofar as they appear random In time) to be obtained with a 
relatively coarse grid. The use of a fine three-dimensional grid, of course, 
requires the use of a large amount of computer time. 

Figure 6 shows a dimensionless velocity component u-| plotted against 
component U 2 (forming a plane In phase space) for one point In physical 
space. Although the behavior here Is much more complicated than that observed 
for the low-order models that are usually used to observe strange attractors 
or strange behavior (here there are not well-defined orbits around fixed crit- 
ical points, possibly because there may be an almost Infinite number of criti- 
cal points), there are similarities. Both the present turbulent results and 
those for the low-order models show trajectories consisting of loops and cusps, 
with frequent changes in the sign of the curvature of the trajectory (e.g., 
Franceschlnl , 1983). (Note that the present turbulent results ultimately 
decay, whereas the low-order models usually do not, since they contain forcing 
terms.) Although there are large changes In the direction of the trajectory, 
particularly In the regions of the cusps, the density of calculated points in 
those regions is very high, so that the numerical results should be reasonably 
accurate. Curves for u^ versus t (the numerical Integrations are with 
respect to t) are. In fact, smooth. Note that u-j and U 2 start out equal 
(on a 45° line) but that their equality Is quickly destroyed when randomness 
sets In. 


^Results from the differenced or spectral forms of the Navler-Stokes 
equations become arbitrarily close to those from the original equations as the 
number of grid points or Fourier components Increases (assuming convergence of 
the numerical method). Theoretically, the Navler-Stokes equations correspond 
to an infinite number of ordinary differential equations or to an Infinite 
number of degrees of freedom. 
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Randomization by strange behavior, or by a loss of hydrodynamic stability, 
almost certainly occurs In the present high-order turbulent results (high order 
here meaning many numerical grid points and thus many equations and degrees of 
freedom). This follows, first, from the fact that such randomness has been 
demonstrated to occur In low-order models (few equations and few degrees of 
freedom), such as that of Lorenz (1963) or Franceschlnl (1983), and second, 
from the fact that with the many more degrees of freedom present here than In 
the low-order models, many more critical points exist, and thus many more 
opportunities for randomization or loss of stability occur. Of course, because 
of the many degrees of freedom here, there will also be randomization by pro- 
liferation of harmonic components (so many harmonic components or eddies pres- 
ent that the Identity of the Individual eddies Is lost and the flow appears 
random). The present large number of degrees of freedom encourages both types 
of randomness, and both very likely occur. 

E. Evolution of mean quantities 

In the results given so far, no correction for discretization error due to 
the finite numerical mesh size was applied. The primary purpose of the present 
work, of course, Is to study the physics of turbulence rather than to obtain 
highly accurate results (possibly unattainable at very high Reynolds numbers). 
For low Reynolds numbers surprisingly good results for the decay can be 
obtained even with coarse grids (fig. 7). At higher Reynolds numbers the 
results, although less accurate, should still be qualitatively correct. Their 
accuracy can be Improved by applying fourth-order extrapolations to zero 
numerical grid spacing (In consistency with the fourth-order numerical differ- 
encing used In the calculations (Delssler, 1981a)). This Is done In lieu of 
subgrid modeling (making an assumption for the eddies smaller than the numer- 
ical grid spacing), e.g., Clark, et al., 1979. The method Is related to sub- 
grid modeling In that It assumes that the subgrid eddies are closely related to 
the calculated eddies but does not require the Introduction of a subgrid eddy 
viscosity (which Is, In effect, a kind of closure assumption). In all of the 
remaining results In this section the fourth-order discretization corrections 
are applied by extrapolating results for three mesh sizes to zero mesh size. 
However, the corrections are negligibly small except at the highest Reynolds 
number. 

1 . Mean-square velocity fluctuations 

Figure 8 shows the calculated evolution of mean-square velocity fluctua- 
tions (spatially averaged) for a series of Initial Reynolds numbers. As the 
Reynolds number Increases (v and Initial length scale xq held constant), 

~2 

the rate of decay of u Increases sharply, as In experimental turbulent flows 
(Delssler, 1979). This can be attributed to the nonlinear excitation of small- 
scale turbulence-like fluctuations at the higher Reynolds numbers. The high 
shear stresses between the small eddies cause a rapid decay. 

2. Microscales and nonlinear transfer of turbulence to smaller eddies 


The development of the small-scale eddies Is seen more clearly In 
figure 9, where the microscale \, normalized by Its Initial value, Is plotted 
against dimensionless time. The microscale Is defined by 
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(51) 


3U^ 3U^ 

For homogeneous turbulence and = 0, X can be calculated from 
equation (15) as 


p Vi 

x 2 cc - M '— 1 (52) 

du^u^/dt 

As the Reynolds number Increases, the small-scale structure becomes finer. The 
microscale decreases until the fluctuation level (Inertial effect) Is low 
enough so that viscous forces prevent a further decrease. After \ decreases 
to a minimum, It begins to grow. (Results for coarser grids were not qualita- 
tively different from these, but the mlnlmums were somewhat higher.) The 
Increase of \ at later times Is due to the selective annihilation of eddies 
by viscosity, the small eddies being the first to decay. Thus, at large times 
only the big eddies remain. It Is this period of Increasing X that Is gen- 
erally observed experimentally In grid-generated turbulence (turbulence 
observed downstream of a grid of wires or bars whose plane Is normal to the 
flow In a wind tunnel). The Increases of X with time observed experimentally 
(Batchelor, 1953, fig. 7.2) are generally of the same order as those In 
figure 9 (doubling the time Increases X by a factor of about 1.5). The early 
period. In which X decreases with time, Is of Interest as Illustrative of 
Inter-wave-number energy transfer. To generate the small-scale structure, 
turbulent energy must be transferred from big eddies to small ones. 

For homogeneous turbulence the equation for the rate of change of 
turbulent kinetic energy (eq. (15)) reduces to 


3_ 

3t 



3U^ 3u^ 


(53) 


That Is, only viscous dissipation contributes to the rate of change of kinetic 
energy, there being no Indication that nonlinear transfer of energy between 
scales of motion Is taking place. There may seem to be a paradox here In view 
of the large transfer of energy to smaller eddies Indicated In figure 9. This 
Is as It should be, however, since energy transfer between wave numbers or 
scales of motion should not contribute to the rate of change of total energy. 
To consider Inter-wave-number energy transfer, two-point equations must be 
used 4 Thus, equation (30) shows that the self-interaction transfer term 
T^ j ( k) In the two-point spectral equation (eq. (26)) has the property that 

T^(k) dK = 0 
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as a spectral transfer term should. The quantity k Is the wave-number vector. 
The spectral transfer term T^(k), or Its Fourier transform 

-3( u^ u j u^ - u^u^ujj/ar^ In equation (19), Is responsible for the generation of 

the small-scale structure In figure 9. Those terms come from the nonlinear 
term -3(u^u^)/3x^ In the unaveraged equation (48). As mentioned earlier, the 

term -3(u^u^)/3x^ produces randomization, as well as spectral energy transfer. 

Although equation (30) shows that T^ can transfer energy between wave 

numbers without contributing to the rate of change of total energy 3u^Uj/3t, 

It says nothing about the direction of the transfer or how Important It Is. 

For that we need calculations such as those In figure 9, which show that sig- 
nificant energy Is transferred to smaller eddies. The energy transfer can 
be thought of as due to a breakup of big eddies Into smaller ones or as a 
stretching of vortex filaments to smaller diameters. In spite of this transfer 
to smaller eddies, experimental results generally show a growth of scale 
(Batchelor, 1953, fig. 7.2). The reason Is that those results are usually for 
the later period shown In figure 9, where, although energy Is transferred to 
smaller eddies, the annihilation of small eddies by viscous action eventually 
wins out. The early period shown In figure 9, and In figure 2 of Taylor and 
Green (1937), Is of particular Interest In that the nonlinear transfer effects 
are truly dominant there; a sharp decrease In scale actually occurs as energy 
Is transferred to smaller eddies. 

3. Dissipation, vortlclty generation, and pressure fluctuations 

The energy dissipation term, the only term contributing to the rate of 
change of kinetic energy for homogeneous turbulence without mean gradients 
(eq. (53)) Is plotted In figure 10. That Is also the mean-square vortlclty 
(eq. (36)), but the two are distinct physical entitles. Although the curve for 
zero Reynolds number, where nonlinear effects are absent, decreases monotonl- 
cally to zero, the curves for higher Reynolds numbers Increase sharply for a 
while and then decrease. Thus, the nonlinear terms In the Navler-Stokes equa- 
tions are very effective vortlclty generators and greatly enhance the dissipa- 
tion at short and moderate times. For long times they appear to have the 
opposite effect, evidently because the turbulence Itself decays rapidly to 
zero. Nonlinear effects, although they do not appear explicitly In the evolu- 


tion equation for u^u^ (eq. (53)), thus alter greatly the evolution by alter- 
ing the dissipation term. 


12 A direct numerical calculation of T^ by Clark et al . (1979) for 
random Initial conditions, and corresponding to the region of Increasing X 
In fig. 8, shows the same thing. Calculated values of T^ are negative at 
small wave numbers (large eddies) and positive at large wave numbers (small 
eddies), so that energy is transferred from big eddies to smaller ones. 
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Figure 11 shows mean-square pressure fluctuations plotted against dimen- 
sionless time. The enhancement of the pressure fluctuations, although not as 
great as that of the vortlclty or dissipation, again Is due to nonlinear 
effects: In this case the nonlinear terms on the right side of the Poisson 

equation for the pressure cause the effect. 

4. Further discussion and summary of the processes In Isotropic turbulence 

Nonlinear velocity and pressure terms do not appear In the evolution 

equation for u.u. (eq. (53)). But we can calculate root-mean-square values 
of the nonlinear terms In the Instantaneous evolution equation (eq. (48)), as 
well as of the linear term. Three measures of the relative Importance of 
Inertial (nonlinear) and viscous effects are shown for a moderate Reynolds 
number In figure 12. The ratio of the nonlinear velocity term to the viscous 
term and the ratio of the pressure to the viscous term In equation (48), 
together with the microscale Reynolds number, are plotted against dimensionless 
time. The terms are space-averaged root-mean-square values. All of those 
measures show a variation from a rather Inertial to a weak fluctuating flow. 

For Instance, R x varies from about 90 to 0.7. This Is a much greater var- 
iation than has been obtained experimentally for a single run. The curves for 
the term ratios lie somewhat below that for R^. They Indicate that except 
at early times the nonlinear Inertial effects associated with velocity and with 
pressure do not differ greatly. 

The Importance of both nonlinear velocity and pressure effects In 
figure 12 Is somewhat paradoxical In view of equation (53), which says that 

neither contributes directly to 3u,u,/3t. The nonlinear velocity effects were 
already discussed In this section; It'was pointed out that such effects should 
not appear In equation (53), since they only distribute energy In wave-number 
space and so do not directly alter the total energy. Although there Is no 
nonlinear velocity term in equation (53), such a term appears In the two-point 

equation for 3u u'/3t. That equation, for the present case. Is obtained from 
equation (19) as 1 
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(54) 


where r Is again the vector extending from the unprimed to the primed point, 
and the pressure terms drop out because of continuity. The last term, where 
the parenthesis Indicates no sum on 1, Is a consequence of the Isotropy of the 

turbulence. The equation for the rate of change of each component of u^u^ Is 

contributed to by the nonlinear velocity term -(3/3r^)(u^u.ju£ - u^u^uj), but 

there Is no contribution from the pressure. The strong effect of pressure 
shown In figure 12 must be contained In higher order equations In the hierarchy 
of averaged equations (moment equations) (Delssler, 1958 and 1960). Thus, 
although two-point averaged equations contain a nonlinear effect of velocity, 
we must consider higher order multipoint equations to obtain an effect of pres- 
sure. Terms In the unaveraged equations shown In figure 12 (averaged over 
space after the solution has been obtained) Include effects of all orders. 
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(Effects contained In the numerical results may, however, be limited by the 
fineness of the numerical grid.) 

Although there is a strong effect of pressure in figure 12, the physical 
significance of that effect Is somewhat elusive In contrast to the effects of 
viscous dissipation and spectral energy transfer. If the turbulence Is aniso- 
tropic, a clear effect of pressure fluctuations Is that they transfer net 
energy among directional components (eqs. (14) and (15) and the discussion 
following those equations). That Is discussed In the following section. If, 

In addition, the turbulence Is Inhomogeneous, pressure can produce a net spa- 
tial diffusion of energy (eq. (15)). Those are evidently the only physical 
effects of pressure fluctuations (at least that we know about). Thus, If the 
turbulence Is homogeneous and Isotropic, as it Is here, It seems reasonable to 
attribute the observed pressure effects In the unaveraged equations to those 
processes. Even though there Is no net Interdl rectlonal transfer or spatial 
diffusion of turbulence when the turbulence Is Isotropic, those processes can 
still be Instantaneously or locally operative. They could, for Instance, cause 
a diffusion of tagged particles. According to figure 12, they have a signifi- 
cant Indirect effect on the evolution of the turbulence. 

From the findings of the present section we conclude that the following 
processes occur In Isotropic turbulence: nonlinear randomization by prolifer- 

ation of harmonic components or by strange behavior, nonlinear spectral trans- 
fer of turbulence among wave numbers or eddy sizes (mainly to smaller eddies), 
spatial diffusion and transfer of turbulence among directional components by 
pressure forces (with zero net diffusion and transfer Into each component), 
generation of vortlclty or swirl, and dissipation of turbulence Into heat by 
viscous action. 13 From this description the life of Isotropic turbulence 
appears Interesting and Includes many aspects. 

V. UNIFORMLY SHEARED FLUCTUATIONS AND TURBULENCE 

In the preceding section the evolution of nonrandom initial fluctuations 
Into Isotropic turbulence was examined numerically. The nonlinear transfer of 
energy to smaller scales of motion, the zero net (but not zero) spatial diffu- 
sion and transfer of energy among directional components, the generation of 
vortlclty or swirl, and viscous dissipation were studied. 

Another Important process Is the production of turbulence by a mean shear. 
Most turbulent flows, both those occurring In nature and those which are man- 
made, are In fact shear flows, where the turbulence Is produced and maintained 
by the shear. Because of the added complexity the nonlinear problem of turbu- 
lent shear flow Is even more difficult than that of Isotropic turbulence. So 
It Is not surprising that little progress has been made In obtaining an ana- 
lytical solution from first principles. An attempt to obtain a numerical 
solution would seem to be In order. 

Conceptually, the simplest turbulent shear flow, although certainly not 
the simplest to produce experimentally (Champagne, et al., 1970), Is one In 
which the turbulence Is uniformly sheared. At least two significant numerical 
studies of that type of turbulence have recently been made (Rogallo, 1981; 


13 Accord1ng to eq. (36), the vortlclty and the dissipation are 
numerically the same, but they are physically distinct. 
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Shaanan, et al . , 1975). In both of those studies random Initial conditions 
with a range of eddy sizes were used. 

In the spirit of the preceding section the present numerical study of 
uniformly sheared turbulence starts with simple determinate Initial conditions 
that possess a single length scale. As In the preceding section, we can In 
this way study how the turbulence develops from nonturbulent Initial condi- 
tions. Again, much higher Reynolds number flows can be calculated with a given 
numerical grid when a single length scale Is Initially present, at least for 
short and moderate times. 

Several Interesting results that could not be obtained In the previous 
work on turbulent shear flow are obtained here. One of the significant find- 
ings Is that the structure of the turbulence produced In the presence of a 
strong shear Is much finer than that produced In Its absence. 

A. Initial and boundary conditions 

In carrying out numerical solutions for uniformly sheared turbulence the 
Instantaneous equations (1) and (3), subject to Initial condition (38) and 
boundary conditions (44) and (47) are used. Since we are considering a 
uniform shear, we let 


dU l 

u i = 6 n dx 2 x 2 


In the Initial condition (38) and 
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In boundary condition (44). For the coefficients In equation (38) we use 
equation (40), where we choose the first set^of signs. Equations (1) and (3) 
are written In terms of the total velocity u . , but we can calculate the 
fluctuating component u^ from equation (4). It should be emphasized that we 
consider here not a sawtooth type of mean velocity profile, but a continuous 
profile In which the mean-velocity gradient is uniform at all points. Even 
with a uniform mean velocity gradient, some local Inhomogeneity Is Introduced 
Into the fluctuations by the periodic boundary conditions. We will not concern 
ourselves with that Inhomogeneity, however, since we can still calculate prod- 
ucts of velocities and pressures averaged over a three-dimensional period. 

Those values are Independent of the position of the boundaries of the cycle. 
Note that for constant uniform mean-velocity gradient and mean pressure the 
last terms in equations (10) and (11) are zero even though the fluctuations may 
be Inhomogeneous (eqs. (12) and (13)). 


The U 2 component of the velocity fluctuation (In the direction of the 
mean-velocity gradient) Is crucial In maintaining the turbulence against the 
dissipation (Oelssler, 1970b and 1972). Therefore, when for brevity only one 
component of the velocity fluctuation Is discussed, that component is chosen 
as U 2 . More will be said about the maintenance of the turbulence later. 
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B. Development of random fluctuations 


Figure 13 shows the evolution of u 2 /\ u q ) at a ^* ed point In space for 

a high Reynolds number, as calculated from the full nonlinear equations. As 
In the preceding section asterisks on quantities Indicate that they have been 
nondlmenslonallzed by using the Initial length scale xg and the 2 k1nemat1c 
viscosity v. Thus t* = (v/Xg)t, x* = x^/Xg, and (dU^/dx^)* = (x^/vldU^/dx,,. 

Again the velocity fluctuations have the appearance of those for a random 
turbulence In spite of the nonrandom Initial condition (eq. (37)). The dashed 

/”2V /2 

curves for u 2 / \ u 0/ are a 9 a ^ n f° r initial conditions perturbed approximately 

0.1 percent. The perturbed curves at first follow the unperturbed ones but 
eventually depart sharply. Although the the curves In figures 13(a) and (b) 
differ considerably In appearance, the perturbed curves In the two figures 
take about the same length of time to break away from the unperturbed ones. A 
very small perturbation of Initial conditions causes a large change In the 
values of U 2 except at small times. On the other hand, the root-mean-square 
values of the velocities change smoothly with time and are unaffected by the 
perturbation of the Initial conditions. These features are characteristic of 
turbulence. (Although the root-mean-square curve In fig. 13(a) appears almost 
horizontal, It eventually goes smoothly to zero when extended.) 


C. Shear-related small-scale structure 


:( u q). 
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tfirttu^ence 


In figure 13 Is the 
(fig. 13(b)) when 
This shear-related 


A striking feature of the curves for 
small scale structure exhibited for sheared 
compared with the structure for no shear (fig. 13(a)). 
small-scale structure Is produced by the term -U|< au^/ax^ In equation (10) 
which, for uniform shear, Is -(dU-|/dx 2 )X 2 au-^/ax-). Equation (10) becomes, 
a constant uniform mean-velocity gradient and a uniform mean pressure, 


for 


2 

3u 1 dU l dU l au 1 a 1 ao 8 U 1 

at = ~ 5 1l dx 2 u 2 ~ dx 2 x 2 ax^ - ax^ ^ u 1 u k^ " p ax^ + v ax^ ax^ 


(57) 


where equation (12) Is used. Note that equation (57) Is obtained without the 
assumption of homogeneity. From the term -(dU-|/dx 2 )X 2 au^/ax-| In 
equation (57), we get the term 


ai) >u »j du a„ u; 

- h; r i "5F7 ■ - r 2 < 58 > 

In the two-point correlation equation (19). For periodic boundary conditions 
x^ dependency Is not present In equations (19) and (58) because averages 
are taken over a three-dimensional period with r^ held constant. This Is 
so even though the periodic boundary conditions may Introduce some local Inho- 
mogeneltles. If we take the Fourier transform of that term, we obtain the 
mean-gradient transfer term T!, In the spectral equation (26). Its effect In 
transferring energy to small-scale components Is similar to that of the non- 
linear transfer term T^ In equation (26) (the Fourier transform of the 
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triple correlation term In eq. (19)). The production of small-scale structure 
by the shear might be thought of as due to stretching of the random vortex 
lines In the turbulence by the mean gradient or to stretching of mean vortex 
lines by the turbulence. 

Although we first discussed a mean-gradient transfer term more than two 
decades ago (Delssler, 1961), the present results give the first graphic 
demonstration of the effectiveness of that term In producing a small-scale 
structure In turbulence. Since that Is a linear effect (when the mean gra- 
dient Is given), we can study It either by the full nonlinear solutions 
already considered In figure 13 (which contain linear as well as nonlinear 
effects) or by linearized solutions. 


Equation (11) becomes, for uniform shear and uniform mean pressure, 


a 2 p 

3X 1 3X 1 ' ~ 3X k 3X k ' 3X 1 ax 2 


(59) 


where equation (13) Is used. As In the case of equation (57), equation (59) 

Is obtained without the assumption of homogeneity. 

D. Some linearized solutions and comparison with nonlinear solutions 

Equations (57) and (59) are linearized by neglecting the terms 
-3(u^U| < )/3X| < and -3 2 (U| < u a )/3X| < 3x^. The numerical solution, with Initial and 
(periodic) boundary conditions given by equations (37), (40), (42), and (43), 
then proceeds as In the nonlinear case. 

We can obtain an analytical solution for unbounded linearized fluctuations 
by using unbounded three-dimensional Fourier transforms (Delssler, 1961). The 
solution does not satisfy constant periodic boundary conditions. Instead of 
working with the averaged equations (Delssler, 1961), It Is Instructive to work 
with the unaveraged ones and use the Initial condition given by equation (37). 
In this case the Fourier transforms must be generalized functions (a series of 
delta functions), but the method of solution Is the same as that In the earlier 
work. Equation (57) for U 2 and equation (59), when linearized, are Inde- 
pendent of U] and U 3 . The solution obtained by using the Initial 
condition (37) Is 

3 

u 2 = u; cos^q n .x - aql^tx^ (60) 

n=l 

3 

p = ^ ^ P n sin ^q n .x - aq^tx 2 ^ (61) 

n=l 

where 
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(62) 
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a = dU.|/dx 2 , q 11 = q^ n + q 2 + q^ , and the a" and q" are given In the 

initial conditions (eqs. (37) and (40)) (with the first set of signs). Mean 
values are obtained by Integrating over all space. For Instance, 
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(64) 


It Is clear from the form of equations (60) and (61) that the solution 
does not satisfy constant periodic boundary conditions. By omitting the term 
-(dU-j/dx2)x23Ui/3X] as well as the nonlinear terms In equations (57) and (59), 
we can, however, obtain a simple analytical pseudo solution that satisfies 
those conditions: 



This solution Is useful for checking the numerical calculations and for study- 
ing the effect of the term (dU-|/dx2)X23U2/3x-| on the fluctuations. 


Velocity fluctuations obtained from linearized solutions (numerical and 
analytical) were plotted (fig. 14). The presence of small-scale structure In 
the curves for (dU-|/dX 2 )* = 4434 and Its absence In those for 
(dU-|/dx 2 )* = 0 are apparent. The curve for no shear (eq. (50)) decays 
monotonlcal ly to zero when extended. This Is In contrast to the nonlinear case 
In figure 13(a) for no shear, where at least larger fluctuations are present. 
The linearized curves for ( dUi /d X 2 ) * = 4434 In figure 14 follow closely 
the nonlinear ones In figure 13(b) for short times. Likewise the linearized 
curves In figure 14 for periodic boundary conditions follow closely those for 
unbounded conditions for short times. For longer times the fluctuations for 
unbounded conditions continue to decay, whereas those for constant periodic 
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boundary conditions grow. Small-scale structure In the curves for unbounded 
conditions Is produced by the term aq^tx 2 In the argument of the cosine In 
equation (60) (a = dUl/dx2). This term arises from the term -ax 2 3U2/3X-] 

In equation (57), as Is evident from Its absence In equation (65), where the 
term -ax 2 3U2/3X] has been neglected. 

For discussing the linearized case for constant periodic boundary condi- 
tions, It Is convenient to convert equations (57) and (59) to a spectral form 
by taking their three-dimensional Fourier transforms. This gives, for u 2 , 
on neglecting nonlinear terms 
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k Is the wave-number vector, and <p_ Is the Fourier transform of u„. Note 
that a finite transform Is used In tne X 2 direction In order to satisfy 
periodic boundary conditions at X 2 /xq = -ir.ir. 

Strictly speaking, equation (67) Is for a sawtooth mean-velocity profile, 
whereas the numerical results are for a uniform mean-velocity gradient. 
Equation (67) should still apply, however, at least for the present discussion 
purposes to points Inside, but not outside, the numerical grid. 

For constant periodic boundary conditions for u^ , small-scale struc- 
ture In the fluctuations or the transfer of energy between wave numbers Is 
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k' In equation (67) 

(eq. 57)) . From Its form 
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that It can produce a complicated inter-wave-number Interaction. 
:y <pl] at each * Interacts with cp, at every other allowable 


produced by the term containing the summation over 
That term Is the Fourier transform of -ax 0 au 0 /ax 
we see 

quantity vpp u u CUV/II iv iiibci uv uj n i wit \f/«p 

K 2 ‘ A difference between the solutions for unbounded conditions and those 
for constant periodic conditions Is that only fluctuations at Integral *2 
are possible when periodic conditions are Imposed, whereas for unbounded 
conditions fluctuations are possible at all values of K 2 - 


The 


Although the linear term -ax 2 3U2/ax-| Is effective In producing 
oscillations, even In the absence of nonlinear effects (fig. 14), the curves 
lack the random appearance of those In figure 13(b). Evidently, as In the case 
of no mean gradients (eq. (50)), the only way we can have a linear turbulent 
solution Is to put the turbulence In the Initial conditions (Delssler, 1961). 
Both -ax 2 au 2 /ax-| and the nonlinear terms In equation (57) are necessary to 
produce the small-scale turbulence In figure 13(b) from nonrandom Initial con- 
ditions. The former acts like a chopper that chops the flow Into small-scale 
components. Although the latter also do that, their most visible effect here 
Is to produce randomization. As In the preceding section the randomization 
might occur as a result of the presence of strange attractors (or, more prop- 
erly, analogous strange behavior) In the flow, by proliferation of eddies or 
harmonic components (with the loss of Identity of the Individual eddies) or by 
both (see section IV for a discussion of these possibilities). 


According to the linearized analytical solution given by equation (60), 
the manufacture of small-scale fluctuations takes place only In the X 2 
direction. Figure 15 shows how this has taken place at a moderate time. 

Figure 16 Is a similar plot for the nonlinear case. The randomizing effect of 
the nonlinear terms Is evident. 


For the nonlinear case, Up'TO/ was a ^ so P 1o ^^ ec * a 9 a ^ ns * x-| (fig. 17). 
The curves show some development' of small-scale structure In the x-| direc- 
tion due to the Interaction of the directional components In the nonlinear 
case. For the linearized flows small-scale structure developed only In the 
X 2 direction. 

E. Evolution of mean quantities with shear 
1 . Cross-correlation coefficients 


7%/2r$/2 

Cross-correlation coefficients u^UjAiJ (uj (1 * j) are plotted 
against dimensionless time for the nonlinear case In figure 17. Although 
~2 ~~2 ~2 

u, = u„ = u at t* = 0, the Initial cross correlations are not zero but are 
all positive and equal. However, because of the apparent randomization of the 

flow u„u„ and u.u„ approach zero as time Increases. On the other hand, the 

CO I o 

values of the turbulent shear stress u.u change from positive to negative 
and remain negative because of the dynamics of the Imposed mean shear. The 
presence of the mean-velocity gradient dU^/dx^ causes u^ to be likely 
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negative when u Is positive, so that u..u , the correlation between the 

two, Is negative: The waviness In the curves In figure 18, as well as that In 
some of the curves In later figures (e.g., fig. 21), Is probably caused by 
nonrandom structure In the flow, possibly that produced by the linear term 
-(dU-|/dx 2 )X 2 au^/ax-) in equation (57) (fig. 14 ). ^ 

2. Growth and anisotropy of the velocity fluctuations 


The evolution of the 


mean-square 
* 

Is plotted In figure 19, where u /lx 2 
adjustment period all of the com 
experiment (Harris, et al. (1977) and 


components of the velocity fluctuations 
2 

(11 c = ( x (/ w ) u m 2. After an initial 

ponents Increase with time, In agreement with 
the numerical results In Rogallo (1981)) 


2 

The numerical results In Shaanan, et al. (1975), on the other hand, show u ? 

~2 

and u 3 decreasing at all times, a difference that remains unexplained. Our 

2 ~2 
u.j component Is the largest of the three, u ? Is the smallest, and u g lies 

~2 

slightly above u 2 , In agreement with experiment Harris, et al. (1977) and 
previous numerical results. 


3. Accuracy of mean and Instantaneous quantities 


2 

The effect of discretization error on the numerical results for u ? Is 

shown In figure 20. Curves are plotted for 16 3 , 24 3 , and 32 3 grid points, 
together with a fourth-order extrapolation to zero grid-point spacing (an Infi- 
nite number of grid points) (Delssler, 1981a, 1981c). The differences between 
the results for 32 3 points and the fourth-order extrapolation are small but 
Increase somewhat at long times. These results appear to Indicate that the 
numerical results given here for averaged values are reasonably accurate. On 
the other hand the three-dimensional spatial resolution Is probably not great 
enough (except at early times) to give accurate spatial variations of unaver - 
aged quantities, other than that they have a random appearance. However, since 
the solutions are hydrodynamlcally unstable, and extremely sensitive to Initial 
conditions, the actual values of the unaveraged quantities are probably not of 
great significance. 

4. Maintenance of the turbulence 


For the case considered In this section (uniform velocity gradient 
dl)-|/dx 2 ) the one-point correlation equation (14) becomes 
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(70) 


where derivatives of averaged values with respect to x-j do not appear 
because averages are taken over a three-dimensional period. This Is so even 
though local Inhomogeneities may occur when periodic boundary conditions are 
used, as discussed earlier. 5 
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Figure 21 shows the evolution of pressure-veloclty-gradlent correlations. 
(Parts of some of the curves are omitted to avoid confusion.) The pressure- 
veloclty-gradlent terms In the one-point correlation equation (eq. (70)), 
together with the production terms, are responsible for maintaining the 
turbulence against the dissipation (given by the last term In eq. (70)). 

~2 ~2 

There are no production terms In the equations for au 2 /at and au 3 /at 

— p ~2 

( «1 1 Uj u 2 aU-j /ax 2 and <s^ u^u 2 au^ /ax 2 are zero). Thus u 2 and u 3 

~~2 

generally receive energy only from the u-j component, whose equation has a 
nonzero production term. Equation (70) shows that to do that, p au^/ax^ + 

p au^/aXj must be positive for 1 = j = 2, 3 and negative for 1 = j = 1. 

Figure 21 shows that Is actually the case for constant periodic boundary 
conditions except for an Initial adjustment period, so that the turbulence Is 

~2 

maintained (fig. 19). The maintenance of the u 2 or u 2 component Is par- 
ticularly critical because If u goes to zero, the Reynolds shear stress 
£ 2 

u,u 2 In the production term of the u. equation (eq. (70)) will go to zero 
and there will be nothing to keep the turbulence going. All of the components 
will then eventually decay. That Is what happens In the linearized analysis 
for unbounded turbulence In figure 21 (see also Delssler, 1961). 

~2 

The nonlinear results for u 2 are compared with various linearized solu- 
tions In figure 22. The same Initial conditions are used for all of the cases 
(eqs. (37) or (38), (40), and (55)). For all of the results, except those for 
the unbounded linearized case (obtained by using unbounded Fourier transforms 

~2 

(eq. (60)), the crucial u„ component eventually Increases so that the tur- 
bulence or fluctuations are maintained. In the unbounded linearized case 

~2 ~2 
u 2 decreases at all times. That was expected, since the u 2 results for 

that case In Delssler (1961, 1970b) (for different Initial conditions) 
decreased at all times. Somewhat unexpected are the linearized results for 
constant periodic boundary conditions, which show that the fluctuations are 
maintained for those cases. Whereas figure 21 shows that In the unbounded 
case the pressure-veloclty-gradlent correlations remove energy from the 

~2 

u 2 component and cause the fluctuations to decay as In Delssler (1961, 1970b) 

the Imposition of constant periodic boundary conditions changes the sign of 

~2 

those correlations and brings energy Into u„, so that the fluctuations are 
maintained. Equation (65), which satisfies periodic boundary conditions, 
shows that, at least when the term ~(dU. ] /dx 2 )x 2 au^/ax^ In equation (57) Is 

— 4 

neglected, u 2 Increases at large times If 2aq!j*q 2 > q n for at least 

one n. 
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5. Spectral transfer terms as stabilizing 

Comparing the linearized case for periodic boundary conditions In 
figure 22 with the corresponding nonlinear case shows that the nonlinear 

~2 

terms have a stabilizing Influence. That Is, the values of u~ Increase 
more slowly for the nonlinear case. Moreover comparing the curve for the 
linearized case with periodic boundary conditions and with the term 
-(dU 1 /dX 2 )x 2 3U2/3X-] In equation (57) missing (eq. (65)) with the correspond- 
ing curve for that term Included shows that the presence of that term also has 
a stabilizing Influence. Since neglect of that term Is equivalent to neglect- 
ing the mean-gradient transfer term T 1 22 spectral equation for 

~~2 

u„ (eq. (26)), we can consider the latter term as stabilizing. Thus both the 
nonlinear spectral transfer term associated with triple correlations T 22 
and the linear mean-gradient transfer term Tp 2 In the spectral equation (26) 

~2 

for u„ are stabilizing. The reason Is that both terms transfer energy to 
small eddies, where It Is dissipated more easily. 


It Is of Interest that the one-point correlation equation for au^u^/at 

(eq. (70)) contains neither a term associated with velocity-gradient transfer 
nor one associated with nonlinear transfer. That Is, both of those processes 

give zero direct contribution to the rate of change of u ^ u j : they only 

change the distribution of energy among the various spectral components or 
eddy sizes. This spectral transfer, of course, still affects the way In which 

u^Uj evolves (fig. 22). Even though equation (70) contains no transfer 

terms, the transfer of energy among the various spectral components of the 

velocity alters the terms that do appear In equation (70), so that 3u u /3t 
Is affected Indirectly. That Is not a small effect! ^ 


The modified linear pseudosolution given by equations (65) and (66) (dash- 
dot-dot curve In fig. 21) Is the simplest solution In which the fluctuations 
can be maintained against the dissipation. In obtaining It the only mean- 
gradient term retained In the equations for U 2 (eqs. (57) and (59), 1=2) 

Is -2(dU-|/dx2)3U2/3X] , a source term In the Poisson equation for the pressure. 
If that term Is also neglected, U 2 decays and, as discussed earlier, all of 
the components of the fluctuations decay. Moreover, as shown In figure 22 and 
already discussed, the term -( dll-| /dX 2 ) X 2 3u^/3x-| In equation (57) Is stabi- 
lizing, so It Is of no help In maintaining the fluctuations. Thus, at least In 
the linearized case, the presence of the source term -2(dU-|/dx2)3U2/3xi In the 
Poisson equation for the pressure Is necessary for maintaining the fluctuations 
That term should play a similar Important role In the maintenance of nonlinear 
turbulence, although in that case It Is hard to separate the linear effects 
from the nonlinear ones. In particular, the role of the nonlinear source term 
In the Poisson equation for the pressure remains unclear, although It may have 
an effect similar to that of the linear source term. 
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F. Return to Isotropy 

Figures 23 and 24 show the approach to Isotropy of nonlinear uniformly 
sheared turbulence when the shear Is suddenly removed. Although the shear 

— 2 

produces considerable anisotropy, the components u^ of the mean-square fluc- 
tuation approach equality upon removal of the shear and remain accurately 
equal. The pressure-veloclty-gradlent correlations In equation (70) are thus 
successful in transferring energy among the various directional components In 

~2 ~2 

such a way that equality of the u. Is produced. We note that u„ continues 
to Increase for a short time afterthe shear Is removed, probably Decause It 

~~2 ~2 

receives energy from both u^ and Ug. 

— 2 

In addition to equality of the u^, zero cross correlations u^u^ (1 * j) 

are required for Isotropy. Figure 24 shows that u^, which Is nonzero when 

the turbulence Is sheared, approaches zero when the shear Is removed, and 
along with the other cross correlations, remains close to zero. The 

destruction of u^u.,, apparently by nonlinear randomization effects, occurs 

over a finite time rather than Instantaneously on removal of the shear. 

Another expected effect of removal of the mean shear Is that the small- 
scale structure produced by the chopping term -(dU-|/dx 2 )x 2 au^/3x-| In equa- 
tion (57) should die out. According to figure 25, that occurs almost Immedi- 
ately when dU-|/dx 2 goes to zero, evidently because of the large fluctuat- 
ing shear stresses between the small-scale eddies. Figure 25 shows, In a 
particularly graphic manner, the effectiveness of the mean-gradient chopping 
term In equation (57) In producing small-scale turbulent structure. 

VI. INHOMOGENEOUS FLUCTUATIONS AND TURBULENCE (DEVELOPING SHEAR LAYER) 

Here, the work Is extended to an Inherently Inhomogeneous developing shear 
layer so that net diffusion, as well as other turbulence processes, can be con- 
sidered. This case Is general enough to Include all of the dynamical processes 
that ordinarily occur In Incompressible turbulence. 

For the Initial conditions we use a three-dimensional cosine velocity 
fluctuation, as before, and a mean-velocity profile with a step. Thus, In 
equation (38) we set 


U Jj = ir6 11 V[sgn(x* - ir) + 1] (71) 

where V Is a constant with the dimensions of a velocity. Equation (71) Is 
plotted against X 2 /xg In the curve for t = 0 In figure 26, where V* = 
Vxq/v , and xg Is again the Initial length scale of the disturbance. For 
the coefficients given by equation (40) we choose the second set of signs. 

With this choice of signs u u does not have to change sign as a result of 
the dynamics of the flow, as'lt did In the last section, and the Initial 
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adjustment period Is eliminated or greatly shortened. If the longer adjustment 
period remained, much of the development of the shear layer would be distorted. 


In carrying out the numerical solution of equations (1) and (3) we use 
boundary conditions (44) and (47), where we let 


(U.) - (U.) 

x*= 2 ir+b* x*=b* 

j j j j 


2 *‘nV 


02 ) 


Equations (1) and (3) are written In terms of the total velocity u. , but we 
can calculate the fluctuating part from equation (4), which, for^the present 
case, Is u^ = u. - $^,U, , where U. Is obtained by averaging u. over x-j 
and x„ for flxea values of x„. The fluctuations are Inhomogeneous In the x ? 
direction, except at t = 0 . 

The calculated evolution of the dimensionless mean velocity U^* = 

(xq/w)U-| ( U 2 and U 3 are zero) Is plotted against X 2 /xq = xj? for a partic- 
ular value of V* = \Ixq/v In figure 26. The results In this section may not 
be as accurate as those In the previous sections because of the presence of the 
discontinuity In the Initial velocity profile, but they should be qualitatively 
correct. The shear layer grows (from essentially zero Initial thickness) 
because of the presence of the turbulent and viscous shear stresses. The ratio 
of turbulent to viscous shear stress (averaged over x-j and X 3 at the cen- 
tral plane X 2 * = ir) Is plotted against dimensionless time In figure 27. 

Except at early times the growth of the shear layer Is almost completely 
dominated by the turbulent shear stress. 


Figure 28 shows the evolution of the Instantaneous velocity component 
U 2 and of the root-mean-square value of U 2 (averaged over the central 
plane X 2 * = ir) . Although the Initial conditions are nonrandom, the evo- 

lution of U 2 has a random appearance, as In the preceding sections. 


On the other hand 


/TV n 

IU 2 j evolves smoothly. 


These characteristics are again 


(~ 2\ W2 

representative of a turbulent flow. The quantity (u 2 J Increases monotonl- 

cally at small times In contrast to the corresponding curve In section V, 
where an Initial adjustment period was present. As mentioned earlier, the 
Initial adjustment period has been eliminated here by using the second set 


of signs In equation (40), so that u u does not have to change sign as a 

n \ 1/2 

result of the dynamics of the turbulence. The decrease In ( u ? I near the 
end of the curve Is caused by a decrease In mean-velocity gradient, and thus 
of turbulence production, at large times (fig. 26). 


As In the case In the preceding section, small-scale fluctuations are gen- 
erated In the Inhomogeneous turbulence In figure 28 by the Interaction of the 
mean velocity with the turbulence. This can be seen by comparing figure 28 
with figures 1 and 13(a), where mean-velocity gradients are absent. One might 
expect this since It has been shown (Delssler, 1981b) that, even for a general 
Inhomogeneous turbulence, a term In the two-point spectral equation for the 
turbulence can transfer energy between scales of motion as a result of the 
presence of mean gradients. 
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A. Inhomogeneous growth of turbulent energy 


A dimensionless plot of turbulent kinetic energy as a function of x * 

2 ^ 

and time Is given In figure 29, where u k u k * = (Xg/v) u k u k . As for all of the 

averaged values In this section, u k u k /2 Is averaged over x^ and x 3 for 
fixed values of X 2 - As time Increases, an Intense concentration or turbu- 
lent energy develops near the plane X 2 /X 0 = •*, where the mean-velocity gra- 
dient Is Initially Infinite. The turbulence Is highly Inhomogeneous. Inhomo- 
geneity, In fact, seems to be the dominant characteristic of the turbulence 
generated In the shear layer. The Indicated Increase In turbulence with time 
Is similar to that obtained experimentally (Brlnlch, et al., 1975). 

B. Turbulence processes In shear layer 

Terms In the one-point correlation equation for the rate of change of the 
turbulent kinetic energy (eq. (15)), which, for the present case, becomes 
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are plotted for t* = 0.000293 In figure 30. As usual, an asterisk on a 
quantity Indicates that It has been nondlmenslonallzed by using x and i>. 

4 3 U 

For Instance (u-|U 2 dU^/dx.,)* = (Xg /v ) u^u^ dU-j/dx,,. The terms that contrib- 
ute most to the rate of change of u^/2 are the production term 

-u.^ dU-j/dx^ the pressure diffusion term (-apu^/a x^)/p and the kinetic energy 

2 2 

diffusion term -(l/2)au k u k u 2 /ax 2 . The viscous diffusion term i>a (u k u k /2)/ax 2 

and the dissipation term -«au /ax au. /ax are small In figure 30. At early 
times, however, when the mean-velocity gradient Is large, the dissipation term 
Is appreciable. 


The production term, whose form shows that turbulent energy Is produced 
by work done on the Reynolds shear stress by the mean-velocity gradient, Is 
largest near the plane xj? = ir, where the velocity gradient Is Initially 
Infinite. The plots of the pressure and kinetic energy diffusion terms show 
that those terms are negative near X 2 = ir and positive away from that 
plane. Thus, they remove turbulent energy from the maximum-energy region and 
deposit It where the energy Is lower. Both diffusion terms therefore tend to 
make the turbulence more homogeneous. 

A comparison of the turbulence diffusion processes with the spectral 
transfer processes and the directional transfer processes arising from the 
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pressure-velocity correlations (sections IV and V) Is Instructive. The spec- 
tral transfer processes remove energy from wave-number (or eddy size) regions, 
where the energy Is high, and deposit It In regions of lower energy. The 
directional transfer processes remove energy from high-energy directional com- 
ponents and deposit It In a directional component (or components) where the 
energy Is lower. The turbulence diffusion processes, as shown here, remove 
energy from regions of space where the energy Is high and deposit It In regions 
of lower energy. The spectral transfer, directional transfer, and turbulence 
diffusion processes tend, respectively, to make the turbulence more uniform In 
wave-number space and more Isotropic and homogeneous In physical space. 

Although one might suppose that turbulence diffusion terms would always 
tend to make the turbulence more homogeneous, that supposition Is not supported 
by all of the experimental data. For Instance, measurements of wall-bounded 
turbulence (Laufer, 1954) Indicate that the pressure diffusion and the kinetic 
energy diffusion terms transfer energy In opposite directions, although the 
total diffusion Is from regions of high to regions of lower energy. On the 
other hand, measurements of turbulence In a free jet (Wygnanskl and Fiedler, 
1969) and In a wake (Townsend, 1949), which are closer to the case considered 
here, seem to support the present findings. 

VII. A STEADY-STATE HOMOGENEOUS TURBULENCE WITH A SPATIALLY PERIODIC 
BODY FORCE 


In all of the cases considered so far the turbulence either ultimately 
died out or Increased In Intensity with time. However, there are many Impor- 
tant cases In which the turbulence, after some time, reaches a statistically 
steady state (e.g., flow In a pipe far from the entrance). Moreover, a discus- 
sion of strange attractors (e.g., Eckmann, 1981, and Ott, 1981) should, 
strictly speaking, be based on a steady-state turbulence; a strange attractor 
Is, roughly, the region of phase space Inhabited by the phase point of a system 
after the Initial transients have died out, where the phase point moves In an 
apparently chaotic fashion. For the decaying turbulence of section IV the 
attractor would then be only a point In phase space. Of course, we could still 
talk about analogous strange behavior, even In an unsteady-state case, as we 
did In section IV-D. 


One way of obtaining a statistically steady-state turbulence Is by adding 
a spatially periodic body-force term (forcing term) F^ to the right side 
of equation (41). A convenient term for that purpose Is 
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(74) 


where the subscript 0 signifies Initial values, the u. are given by 
equations (38) and (40) with U^ = 0 and p by equation (3), and c Is a 

constant. The first set of signs Is used In equation (40). Equation (74), 
which Is time Independent^ Is used for F^ at all times. For c = 1, the 
the quantities u^ and p, as calculated from equation (1) (with F^ added 

to the right side) and equation (3), do not change from their Initial equation 
(3), do not change from their Initial values. To Introduce some Initial time 
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dependence, we set c = 1.05. The boundary conditions are taken to be 
periodic, as In section IV. 

Calculated results for„th1s case are plotted In figures 31 and 32, where 
t* Is again equal to {v/xq )t and xq Is the Initial length scale. 

Figure 31 shows the time evolution of u^ and ( u.j I at a point away from 
the center of the numerical grid, where as before, overbars Indicate space 

r2\ 1/2 

averages. The values are normalized by dividing them by (u n I , the Initial 
/^\ 1/2 /^\ 1/2 [~ 2 \^ /2 

value of I u, 1 = (u ? ) = (u 3 J . Since we are Interested In steady-state 

solutions'll Hong times It Is'necessary, to obtain reasonably accurate results, 
to use a lower Reynolds number than In the preceding cases, where shorter-time 
transient flows were considered. 


Figure 30 shows that, for 0 < t* < 0.17, the flow Is essentially laminar 

18 1 

1/2 


with small fluctuations of u^ . Then for 0.17 < t* < 0.18 there Is a rather 


sharp transition from laminar to turbulent flow, as 


(*?) 


Increases. For 


t* > 0.18 the turbulence Is statistically steady state, as Indicated by the 

/ToO /2 


nearly constant value of 


(=:) 


Curves for u 2 and u 3 are similar to 


those for u, , Including the same location of the transition region and nearly 

/^ 2\ 1/2 /^ 2\ 1/2 / 3 \ 1 /2 

the same values for (u^j and (u^J as for (u^j 


After the transition region (t* > 0.18) the flow appears to lie on a 
strange attractor, since it has the following characteristics: 14 first, a 
volume In the phase space of our system decreases with time, since the Navler- 
Stokes equations describe a dissipative system and phase-volume shrinkage can 
be shown to occur for the Navler-Stokes equations. 14 This Implies that an 
attractor exists for our system. Second, the chaotic appearance of the veloc- 
ity components (figs. 31 and 32) Indicates that the attractor Is strange. 
Finally, the fact that transients have died out for t* > 0.18, leaving a sta- 
tistically steady state (fig. 31), Indicates that beyond the transition region 
the phase point Is on the strange attractor. 

Figure 32 shows the projection on the u„-u 3 plane (at the center of the 
numerical grid) of the trajectory of the phase point as It moves on the strange 
attractor. As In figure 6 the trajectory consists of loops and cusps with fre- 
quent changes In the sign of the curvature, but unlike figure 6 the trajectory 
does not, of course, tend toward a point. The cusps might be considered as 
loops with very small or zero radii. Also, as In figure 6, randomization Is 
very likely associated with the large number of harmonics (eddy sizes) present, 
as well as with the strange attractor or strange behavior. 


14 These characteristics, as well as the possibility of obtaining a 
steady-state turbulence with periodic boundary conditions by slightly modifying 
the existing program, were pointed out to the author by R.J. Delssler. 
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VIII. CONCLUSIONS 


From the present review It Is concluded that the nonlinear and linear 
processes In turbulence can be profitably studied numerically. The results 
show that, at least at higher Reynolds numbers, an apparently random turbulence 
can develop from nonrandom Initial conditions. The numerically calculated 
turbulence Is not numerical hash, since a large number of time steps correspond 
to each fluctuation. For both sheared and unsheared fluctuations the structure 
of the Navler-Stokes equations Is such that turbulence can develop even when 
the/lnltlal flow Is nonturbulent. This Is Indicated by the appearance of the 
Instantaneous velocity fluctuations and by the sensitivity of those fluctua- 
tions (and the Insensitivity of average values) to small perturbations In the 
Instantaneous Initial conditions. The randomness appears to Increase as the 
numerical mesh size decreases. Moreover, the two-time velocity correlation 
becomes small as the time between the occurrence of the two velocities 
Increases. In addition, for no mean shear the correlation between any two com- 
ponents of the velocity becomes small as the time Increases, as a result of the 
randomization. This correlation Is not small Initially, even though the three 
components of the mean-square velocity fluctuation are equal at early as well 
as at late times for the Initial conditions chosen. Also, calculated velocity- 
derivative skewness factors for no mean shear appear to be of reasonable mag- 
nitude when compared with those for Isotropic turbulence. Thus, except In the 
Initial period the results for no mean shear evidently give a reasonably good 
approximation to Isotropic turbulence. 

The source of the observed randomness may lie In the presence of strange 
attractors or, more properly, of analogous strange behavior (Monin, 1978) In 
the phase space of the system, as well as In the occurrence of a large number 
of eddies or harmonic components (large number of degrees of freedom). It 
appears that no conclusions can be drawn as to the relative Importance of the 
two processes, but both probably occur. (A strange attractor Is a region In 
the phase space of the system to which solutions are attracted and In which the 
phase point moves In an apparently chaotic fashion. It can occur even with a 
small number of degrees of freedom.) Roundoff errors appear not to be a sus- 
taining cause of the randomness; a large decrease In roundoff errors did not 
appreciably affect the turbulence level or the randomness of the fluctuations, 
although the Instantaneous values were different. Thus, the affect of a large 
decrease In roundoff errors Is similar to that of a small perturbation of the 
Initial conditions. Roundoff errors may In some cases affect the transition 
to turbulence. The present turbulent solutions bear some similarity to those 
for low-order models In that both have trajectories In phase space that consist 
of loops and cusps, with frequent changes In the sign of the curvature of the 
trajectory (fig. 6). Moreover, with the results from the low-order models In 
which apparent randomness appears with as few as three degrees of freedom 
(e.g., In the Lorenz equations), the turbulence observed to be manufactured by, 
the Navler-Stokes equations should perhaps not come as a surprise. 

At early times the calculated nonlinear transfer of energy from big eddies 
to small ones Is almost completely dominant and causes a sharp decrease In the 
size of the microscale. This has not been generally observed experimentally 
or analytically because the period usually studied Is for later times, where 
the annihilation of small eddies by viscous action causes the scale to grow, 
even though energy Is being transferred to smaller eddies. This later period 
of scale growth Is also observed In the present results. 
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The nonlinear terms In the equations of motion, besides transferring 
energy among eddy sizes and producing randomization, are very effective vor- 
tlclty generators and Increase the dissipation and the rate of decay. The 
Increased rate of decay Is a result of the nonlinear transfer of energy to 
smaller eddies; the small eddies decay faster than the big ones because of the 
higher shear stresses between the small eddies. Calculation of (averaged) 
terms from unaveraged equations of motion shows, as might be expected, that the 
flow Is dominated by nonlinear Inertial effects at early times and by viscous 
effects at later times (fig. 12). The nonlinear effects are associated with 
both velocity and pressure terms In the unaveraged equations of motion, even 
for Isotropic turbulence. Since the one- and two-point averaged or correlation 
equations for Isotropic turbulence do not contain pressure terms, the effects 
of pressure observed for the unaveraged equations must be contained In higher 
order averaged equations. The Infinite hierarchy of averaged equations should 
contain all effects, as do the unaveraged equations. The only physical proc- 
esses associated with pressure (that we know about) are Interdlrectlonal 
transfer and spatial diffusion of turbulence (eq. ( 14 )). it thus seems rea- 
sonable to attribute the observed pressure effects In the unaveraged equations 
to those processes. Even though there Is no net Interdlrectlonal transfer or 
spatial diffusion In Isotropic turbulence, those processes can still be locally 
operative. 

The processes occurring In Isotropic turbulence thus Include the follow- 
ing: nonlinear randomization, nonlinear spectral transfer (mainly to smaller 

scales of motion), zero net (but not zero) spatial diffusion and transfer of 
turbulence among directional components, generation of vortlclty or swirl, and 
viscous dissipation. 

If a uniform shear Is present In the flow, we have, In addition to these 
processes, production of turbulence by the mean-velocity gradient, net transfer 
of turbulence among directional components by pressure forces, and linear 
spectral transfer among scales of motion by the mean gradient. The last of 
these processes results In the production of small-scale fluctuations In the 
flow. This can be attributed to a mean-gradient transfer term In the spectral 
equation for the velocity fluctuations (eq. (26)). Although we first discussed 
that term over two decades ago, the recent numerical results considered herein 
give the first graphic demonstration of the effectiveness of that term In gen- 
erating a small-scale structure In the turbulence. However, the small-scale 
fluctuations produced by that term alone (linear solution) are essentially 
nonrandom. Evidently, the only way we can have a turbulent linear solution, 
either with or without mean gradients, Is to put the turbulence In the Initial 
conditions. To produce the small-scale turbulence from nonrandom Initial con- 
ditions observed herein for shear flow, the presence of both the linear mean- 
gradient transfer term and the nonlinear terms In the equations Is necessary. 
The former term, or Its equivalent In the unaveraged equation (57), acts like 
a chopper that chops the flow Into small-scale components. The latter terms, 
while they also produce small-scale components, act most visibly here as 
randoml zers . 

In all of the uniform-shear cases calculated with constant periodic bound- 
ary conditions, Including both linear and nonlinear flows, the pressure- 
veloclty-gradlent correlations are successful In distributing energy among the 
directional components, so that the turbulence or the fluctuations are main- 
tained. This Is In spite of the presence of a production term In the equation 
for only one of the components. Both the linear mean-gradient transfer term 
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and the nonlinear terms mentioned In the preceding paragraph have a stabilizing 
effect. That Is, they cause the fluctuations to Increase at a slower rate. 

The reason Is that both terms transfer energy to small eddies, where It Is 
dissipated more easily. It Is shown that, at least for the linearized solution 
with constant periodic boundary conditions, a mean-gradient source term In the 
Poisson equation for the pressure Is necessary for maintaining the fluctuations 
against the dissipation. That term should play a similar Important role In the 
maintenance of nonlinear turbulence, although In that case It Is hard to sepa- 
rate the linear effects from the nonlinear ones. For the linearized unbounded 
solution (obtained by using unbounded Fourier transforms) the fluctuations 
decay, as expected from earlier results. 

When the mean-velocity gradient Is suddenly removed, the turbulent shear 
stress goes to zero In a finite time, and the velocity-pressure-gradient cor- 
relations cause the turbulence to attain the Isotropic state. The Intensities 
of the directional components become and remain equal. In addition, the small- 
scale structure produced by the mean-gradient transfer term quickly vanishes 
(fig. 25). Figure 25 shows, In a particularly graphic manner, the effective- 
ness of the mean-gradient chopping term In equation (57) In producing small- 
scale turbulent structure. 

For a developing shear layer the turbulence Is Inhomogeneous and, In 
addition to the processes considered so far, a net spatial diffusion of turbu- 
lence occurs. The thickness of the shear layer, which Is Initially zero, 
Increases with time because of the presence of turbulent and viscous shear 
stresses. Except at very early times the growth of the shear layer Is almost 
completely dominated by the turbulent shear stress. As time Increases, an 
Intense concentration of turbulent energy develops near the plane where the 
mean-velocity gradient Is Initially Infinite. The turbulence Is highly Inhomo- 
geneous. The calculated turbulence production Is always positive, and Is 
largest near the plane where the velocity gradient Is Initially Infinite. The 
pressure and the kinetic energy diffusion are negative near that plane and 
positive away from It. Thus, they remove turbulent energy from the high-energy 
region and deposit It where the energy Is lower. Both diffusion processes 
therefore tend to make the turbulence more homogeneous. 

A comparison of the various transfer and diffusion processes occurring In 
turbulence Is of Interest. The spectral transfer processes remove energy from 
wave-number (or eddy size) regions where the energy Is high and deposit It In 
regions of lower energy. The directional transfer processes remove energy from 
high-energy directional components and deposit It In a directional component 
(or components) where the energy Is lower. The turbulence diffusion processes 
remove energy from regions of space where the energy Is high and deposit It In 
regions of lower energy. The spectral transfer, directional transfer, and 
turbulence diffusion processes tend, respectively, to make the turbulence more 
uniform In wave-number space and more Isotropic and homogeneous In physical 
space. 

By adding a spatially periodic body-force term to the Navler-Stokes equa- 
tions, a solution Is obtained In which the flow first passes through laminar 
and transltlon-to-turbulence stages. The turbulence then quickly settles down 
to a statistically steady state. In this last stage the flow appears to have 
characteristics corresponding to those of a strange attractor. 
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Initial conditions use eq. (40) with 

first set of signs 



t* 


(a) x* = x* = 9ir/8, x* = 3ir/8, for unaveraged fluctuations. 

(b) x* = x* = x* = ir, for unaveraged fluctuations. 

I C. J 

Figure 1. - Calculated evolution of turbulent velocity fluctuations (normalized by initial condition) 

fiV /2 

for a high Reynolds number ^u Q j x Q /v = 2217. No mean shear; root-mean-square fluctuations 
spatially averaged; 32 3 grid points. 
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(a) 8 3 grid points. 

(b) 16 3 grid points. 


(c) 32 3 grid points. 


Figure 2. - Effect of numerical mesh size on calculated evolution of velocity fluctuations. No mean 


~ 2\ ]/2 


shear;! u^ x /\> = 2217; x* = ir (at grid center) 




5 



*#_ f * 

' *0 

(a) tQ = 0. 

(b) t 0 = 0.00813. 


Figure 3. - Calculated correlation coefficient for velocities at dimensionless times to and to 



t* 

Figure 4. - Calculated correlation coefficient for two velocity components plotted against dimension- 

, , . 

less time. No mean shear^UQ J xq/v = 2217; 32 3 grid points. 
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t* 


Figure 5. - Calculated evolution of velocity-derivative skewness factor. No mean shear; 



Figure 6. - Calculated trajectory of phase point projected on u 1 ~ u 2 P 1ane x* = x* = 9ir/8, 

x* = 3ir/8, and 0.00236 < t* < 0.0108. Arrows indicate direction of time. No mean shear; 

3 

3 

32 grid points. 
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Figure 7. - Effect of numerical mesh size on evolution of u at low and moderate Reynolds numbers. 

u u ~2 ~2 ~2 ~2 

No mean shear, u = u^ = u 2 = u^. 



Figure 8. - Calculated evolution of mean-square-velocity fluctuations (normalized by initial value) 

2 2 2 2 

for various initial Reynolds numbers. No mean shear; u = u^ = u 2 = iij ; extrapolated to mesh 
size. 
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Figure 9. - Calculated evolution of microscale of velocity fluctuations (normalized by initial 
value) for various initial Reynolds numbers. No mean shear; extrapolated to zero mesh size. 



2 

Figure 10. - Calculated development of mean-square-vorticity fluctuations u or dissipation c 
(normalized by initial value) for various initial Reynolds numbers. No mean shear; extrapolated 
to zero mesh siie. 
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Figure 11. - Calculated evolution of mean-square pressure fluctuation (normalized by initial value) 
for various initial Reynolds numbers. No mean shear. Extrapolated to zero mesh size. 



Figure 12. - Three measures of relative importance of inertial and viscous effects plotted against 

“2 1/2 

dimensionless time. No mean shear; u^ x^/v = 69.3; i = 1, 2, or 3; extrapolated to zero mesh 
si ze . 
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Figure 13. - Effect of uniform shear on calculated evolution of nonlinear turbulent velocity 

fluctuations (normalized by initial value) for a high Reynolds number lu^ I x^/v = 1108 

Root-mean-square fluctuations are spatially averaged; 32^grid points; x* = x* = 9v/8, 
x* = 3ir/8 for unaveraged fluctuations. c 
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Figure 14. - Calculated evolution of linearized velocity fluctuations (normalized by initial value). 

H\ yZ 3 

Uq ) x Q /v = 1108; x* = x* = 9ir/8; x* = 3ir/8; 32 grid points. 
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Figure 16. - Nonlinear solution for u^u^ j plotted against x*. x* = 9ir/8; x* = 3ir/8; 
(dU 1 /dx 2 )* = 4434;^) x q /v = 1108; 32 s grid points. 
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Number of 
grid points 
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Figure 21. - Calculated evolution of pressure-velocity-gradient correlations. (dlt^A^)* = 4434; 

( 2\ 1/2 3 

I Uq 1 Xq/ v = 1108; 32 grid points. 
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Figure 22. - Evolution of u„ 
(~ 2\ U2 

Uq x 0 /v=H08. 


for various linear and nonlinear solutions. (dU-|/dx 2 )* = 4434; 



Figure 23. - Calculated approach to isotropy of uniformly sheared turbulence upon sudden removal of 

ni /2 3 

the shear, lu* ) x Q /v = 1108; 32 grid points. 
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0 .2 .4 .6 .8 tor 1.2 1.4 1.6 1.8 2 n 

x 2 " x 2^ x 0 

Figure 26. - Calculated development of shear layer mean-velocity profile with dimensionless time. 

rr \ 1/2 3 

(u Q ) x Q /v = 554; V* = 2216 in eq. (71); 32 grid points. 
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Figure 27. - Calculated time variation of ratio of turbulent-to-viscous shear stress for developing 
shear layer at x$ = ir. ^ J »q/v = 554; V* = 2216 in eq. (71); 32 3 grid points. 



0 2 4 6 8xl0 4 


Figure 28. - Calculated evolution of turbulent velocity fluctuations (normalized by initial value) 
for developing shear layer. Unaveraged fluctuations are calculated at center of numerical grid 

(x* = ir). Root-mean-square fluctuations are averaged over x* and x* at central plane 
1 — wo 13 

* ( 2 V 2 

X 2 = ir. J fy/v = 554; V* = 2216 in eq. (71); 32 3 grid points. 
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x 2 " V x 0 

Figure 29. - Development of dimensionless kinetic energy profile with dimensionless time for develop- 
(~2\ yi 

ing shear layer. ^uqJ xq/v = 554; V* = 2216 in eq. (71); 32 3 grid points. 
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Figure 30. - Plot of terms in one-point correlation equation for kinetic energy (eq. (73)) for devel- 

rz\' /z 

oping shear layer. \uq ) xq/v = 554; V* = 2216 in eq. (71); t* = 0.000293; 32 3 grid points. 
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Figure 31. - Calculated evolution of turbulent velocity fluctuations (normalized by initial condi- 

tion) with a spatially periodic body force. I u Q I x Q /v = 138.6; x* = x* = 9ir/8; = 3ir/8 for 

unaveraged fluctuations. Root-mean-square fluctuations spatially averaged; 32^ grid points. 



Figure 32. - Calculated trajectory of phase point projected on u^-u^ plane at numerical grid 

center with a spatially periodic body force. 0.26 < t* < 0.32; 
points. Arrows indicate direction of time. 


-i \ 1/2 

u 0 ) x Q /v = 138.6; 323 grid 
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